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ABSTRACT 


The problems of tension, simple shear and torsion are 
solved using a hypoelastic constitutive equation similar to the 
Prandtl-Reuss equation for an elastic-plastic material. It is shown 
that the hypoelastic solutions give underestimates for the force vari- 
ables of the classical elastic-plastic solutions in these particular 
problems. 

Simple extension and torsion of a circular and Square bar 
are investigated experimentally for commercially pure aluminum. A 
work-hardening function is approximated from the results of the tension 
test. The torsion problem is then solved for both the circular and 
Square cross section using hypoelastic theory. The theoretical and 


experimental results are compared for the work-hardening material. 













TIAA TAA 





ans notetot bad vaste sfamte ,notans? 76 ematdorvg SAT 
add ot veftmbe notssups: svtsus (p2ne5 sttasleoqyd & pneu bavfoe 
nwode ef +! .tetrsten atten g-strasts: 16 yot notteupe 2euah-[tbnsv9 
-hv6v sanot silt 10% zetsmbtesvobru ovtp anotiwioe ottesisoqyd ems tana 
‘eluoriveg g2ad nt enofsuloz otteslq-ottesle [eatzests sit to 2afds 
ame dem 

isd oveupe brs ‘eiusit> 6 to nofte1ot tne notenetxs sfqute 
A .muntmwls sw yiistoremmoo sot yl Tatneniyeqxe betepiieevat ste) 
fotenst ssi3 to ativeer siz mor? betsmixovdgs 2? nettonu? patnebrad-! 
bns vsfuarts sag ritod yo? bavioz ned ef mafdovg notevad SAT  tzeF 
bine [sattevoord ant ,yioont attest sequel prfev notioe2 22049 syseup2 — 


Setvstem pninebysi~avow end 107 boysqmoo s¥6 2dfuesy (sdnemPyaqxs 





ACKNOWLEDGEMENTS 


The author wishes to express his gratitude to Dr. J.B. 
Haddow for the guidance and supervision of this thesis. 

Appreciation is extended to the National Research Council 
and the University of Alberta for the financial support received 
during the course of this research. 

The experimental test specimens were prepared by members of 
the Mechanical Engineering shop. Mr. R. Haswell assisted in all stages 
of the experimental work and Mr. F. Christopher was of great value 
during the instrumentations. Typing of the thesis was done by Miss 
H. Wozniuk. 

The author would also like to thank his wife, Sigrid, for 


her patience and encouragement. 


iv 














2TH IMIG I INONADA 


O.t .10 of sbuditerp et 2esrqxe of eonetw vOnsiE ont 
.2tagiit efit to notetvraque bas sonsbtup ent 107 wobbelt 

[fanved dove9298 [enofssi ont of babnsdxe et notterosyggA |S 
bavfene tvoqque Tstonsnt? ont v0? s3NsdfA to Yiteyevialoedd Bas | 
.foves297 2tdt to senwoo ont enkvwb 

+0 2vedmem vd boteqe1q 919W ensmiaaqe test fstnsmiyegxe snT od 
eopate ffs nt betzrezs lowes .A wiM .qode patvesnipnd Testnedos sag ae 
sulav teoxp to 2ew vatqotertydd .7 . 1M bas diow lednomimisaxe ait 4 
22tM yd shob 2sw 2tesnds srt to onfaut .enotisinamintent ont ¢ 

| a 


yet ,.bhapt2 ,stiw 2rd snes ot vif o2fs bi wow “oftus oAT 1 Ae 


.Inemepsyueons brs sonsiteq, 19d 


S 


CHAPTER I 
CHAPTER II 


CHAPTER III 


TABLE OF CONTENTS 


LNT ROTIDETIEON FARR Gate, cre sstesetarstesetatateteteratate ete tu tefatatetatavete 


THE RELATION BETWEEN HYPOELASTICITY AND 
ROE Bier Syiate emer te, deste <a wlnch “e & “dinrs “aca Wc es 


2.1 The General Hypoelastic Constitutive 
EQUAL TON) © ere ete AIRS ARSE RING wg cas 


2.2 An Elastic-Perfectly Plastic Constitutive 
Equation GACERAAELON..Q7, G16 Oy nanone cen... 


2.3 A Hypoelastic Constitutive Equation Similar 
to the Prandtl-Reuss Equation ........... 


2.4 A Hypoelastic Constitutive Equation Suit- 
able for Some Work-Hardening Materials .. 


SOLUTIONS OF SOME HYPOELASTIC PROBLEMS AND COM- 
PARISON WITH CORRESPONDING ELASTIC-PERFECTLY 
PEAS PE GESOL UL LONGMIRE CS ies 6 ak oie wie Gis ohdetene a «avers 
Om le AINE TE XUCUSMOM ear isi, cee sae 0 ce sok ao chin 
CRPer AS MTN CPONCAT SHER Os Riots Sane SST SEAS EES 
3.5, MNOVSTONIG Tees. be eC emergeni ag... 
oeo7ds wotatement.or-ctne Problem....¢..<.. 
JeoeDe ClaSSiGadimelaStic LOYSION ..:sa6.: 


acorn HyPOGlasS tier (OVS1 ON. 6.6 nc cc ance ae 


3.3.d Numerical Procedure for a Square 
GrosS=Seetionons bib. YS. PaveS ee. 


3.3.e Results for Torsion of a Square 
BAY’: wenusketegs despises «5 cece eee eee 


3.3.f Hypoelastic Torsion of a Circular 
DOVE 6 utctese, tale) cers cael iene eae aaa eats 


10 


17 


22 


24 
24 
26 
36 
36 
36 
40 


5] 


Da 


60 









ad ' 
Te THI 
¢ ‘ r + be Ti 7 
» i } 3G - ‘4 { | Biel aA 
‘ yT 
cciee so ot e's wt Fs e ’ I if i { 2A J9 
if b i : [ $ 
¢ om 
w Agr 4 
} 4 =m ; 1 si 4 ss S 
e 


~» NOrs5ups 


2 ‘ aus " fi * ££ 
ae ) i j J 4&4 %d +3 26 i500 : a) Gu oa 
f ‘ P - fe - = - 
‘ 0 H-)IDAbYT Set of 


¢ 
~~ 
= 

> 
be 
ca 
fy 
_— 
m4 
_ 
wa 


4 
» 
i= 
< 
> 


4 
> t... 
| a Cc 
od utd 


ee 
™~ 


oy 
= 
a) 
u 
2 
= 
~w 
i 
| ep 


ot 7 ’ 4 am > * rroivk . ~ c- 
-e+ coca H Ct j U jnoemetate Dee 


Ve . * r? ae fest ae fy a , 
np Pe noreyol prtesla. Isoteesid dst.t 





oh 4 Al i ‘ 
. Qe ee ee mofeyoy attesf[ooqvH aib.b- - 7 














_ St6upe 5 VOT stubsD079 Tsotnotun b. ~ £ 
We OVE bad an dhin ene ay notdoee+zeny) |p 





CHAPTER IV 


CHAPTER V 


TABLE OF CONTENTS (continued) 


EXPERIMENTAL CONSIDERATION OF A WORK- 
HARDCNINGs MATER LAL et ettete cn. enlu se ea ee ncay oe ee 


eee OOS ON ALCS US oie eres, cis bie das cs wie od lev ws 


4, 
4. 


4. 
4. 


] 


Te 


] 
] 


al 


2 


4 
re) 


iG 


Experimental Apparatus ........... 
ThGsiese SPECIMENS: canwss ce owen ss 


Calibration and Experimental 
GOR CUULGMIT ES «ss choosy eeironence 


Calibration of the Dynamometer.. 


Calibration of Twist Measure- 
(NGI Rete sk esac seis cs cece wa 


Experimental Procedure ......... 
Results of the Torsion Test ...... 


Coma aOmmteS lect c.f cicin ais soot 4 «os 


ARO MLOWS LOM LCS Cemeeeeiett Sy secs oo ccc este vice shee wala 


SOLUTION OF SOME WORK-HARDENING PROBLEMS AND 
COMPARISON TO EXPERIMENTAL RESULTS ........... 


5.1 Determination of the Work-Hardening 
BUG CMON: cere etk aris cee sien ek see ales 


Bc MEA) LOS Caterer et. 2 cic. cs <5 6 o6 4 sb ecdy oe 


POLS LON, OL Oe G mr CULO aD AG. ot sas ecco ss cet 


Beas eLOGS ON VOT pamsUUd CORD Co oc c/s tro ieecs,stegers aretans 


5.5 An Elastic-Plastic Constitutive Equation 


for a Work-Hardening Material 


vi 


114 


114 
117 
1g 
120 


| ae 














7 my : ' 
- ig 
ih ae “~ i 7” 93) Poll Be Fil Vor ;’ 10_ 230A 


-N8OW A TO HOTTANIOLSHOD. JATMAMIAIKS ==“ VT AITIAHD 
a eccevveccssess JATRSTAM OMIMSCGHAH 


see eevee tees 


Jn wean Cee ‘257 norevoT: 7,8 : 
ary B: + | ' 3 UA. [. me 
7 ig2 325T SdT ~ 8.7.8 
DNS NObIETGI [60 t.1,8 
( a Ar - teoseee S1UDSDOTVY 
_ 4 = 
~ iO i TUT 16) F ct wh 
Mm FOTW + ( J > i 6.) d £ : de 
Se ae italien ae an ee J ; 
grub dnsurysaxd 3.8.7.8 
es! NOrevol sit to erivesh #.f.4 
PO Re as ee [ ry ). 2 h 
ee, ee ae iesi nofenel §.4 
if oP s | : AV iad HAW = OW IMO: 40) VvOTTUJOZ - y 7 TIAr ; 
iT nH1ISIG iA ro ‘) ‘ ar AT ~ a> eae a 
JU238 JATVAMIFS9XS OT MOZTHAIMGD s 
DM Fis 18 ke Oe | O it) a4 3.5 90 {2 
~ = 
75 = —~— 90 @ eee, oe 6 eels eee « . . . ’ - fit 2 Th s 
ert ; : 
GT 5 1c rege ee ise Wed afqmt2 S.2 
7 a 
4 Tie, -eereeees “6G YSINOVTD 6 tO MOTe10T £.2 [ 
; 4 
' QS -* -~e* eee eo eweneeae NBG Sisuoe -6 bie) nor240) 4.2 : 










notteups svigusittieno) stseeltatses! aA 2.2 
Sean oo f AN mos a tna tcH~AYOW 6 vor . 
7 a | 





i] 


CHAPTER VI 


BIBLIOGRAPHY 


TABLE OF CONTENTS (continued) 


CONCLUSION ©) scares 5 ov nk ts ecsieialeia eels ove 
Poet IMAL TONBOREERTONS onto cei ees a sles onus 
6.1.1 Errors in Numerical Procedure .... 
Deal eee EXDebINieica | EYYOYS! «xisa 05 os 2 6 sc 5 
ee) SEUSS OMSEEETOLS go srtly o otey in se als & ee 


Pe Tee TSCUS SI OMMGUERGSULUS Coc bac dc cece es bce 


eoeerer eee eee eee eee eee eee oes eoee eee eevee oeee ee eeeeee 


vii 


Page 
134 


134 
134 
14] 
142 
143 


eye 


; 
“+eeeee re eee eee ane® ‘tie * ere oa 8 o> VOT2UIINOD 


ey Ar. e103 49 or kamite3 

. sxubso079 [sof remun ni avov"a r.4.a 
Jaen eros14 fesnsnirsgxs $.j/.0 

doe Ate ee 2v0v1d to nofeevserd 


Pe nth tee + ed fU 259" to nofzevo2td 


a 









$.8 
£.3 


eee eeweeeeeew te eevee eee + **-* ee #r ee f# @ m8 ® 





> 


> 


FP FSP FSP LP FP LP LFS +f 


o13 
14 


LIST OF TABLES 


Results for Elastic-Plastic Torsion of a Square 
TENDS save Bree eeu SoS. by os RR Re 


PO ROUCE Go Lt 0 boven Ol eee oo ccs occa sine oh dang aie suas 


Calibration of Angle of Twist for Circular 
SiOx petit) ae St Ae | Se eee es 


Calibration of Angle of Twist for Square Specimen. 
Experimental Data for Torsion of Circular Bar .... 
Results/ for Circular Cross-Section: ciscsscerissis: 
Experimental Data for Torsion of Square Bar ...... 
Results for Square Cross-Section ............0000. 
sunmaey Of #Resuits! for Circular: Bar Veiee eo 390. e 
Summary of Results for Square Bar ................ 


Calibration of Angle of Twist for Comparison 
SHECIMEN 9 ied ct aweeneet digs estieter resis es 


Experimental Data for Comparison Test ............ 


Calculation of Angle of Twist from Spider 
WOUCAGE TREACINGS OU. Fimo we lee eer. PET. SUL RS ere. 


Experimental Data for Tensile Test ............... 
Summary of Results for Tensile Test ............6- 
Estimated Slopes of Tensile Results ............6. 
Resutts “for STors ion ‘oreGircutar -Bar +e. centers: 
Estimation of Maximum Errors in Slopes ........... 


Estimation of Maximum Errors in Moments .......... 


viii 


100 
101 
127 
128 
149 
150 


———-— 





















yyeupe ‘io norevoT zsttesld-atteatls yor etiwveah 
eee 6 6 4 eget a oerveoeweeve eee ea eeee "5a 


avons oe worterdT 162 supyoT 


ely io} t2iwT to ofond to nofPsesdhlso 
*@seteetie*e va <ses a» NOmMtosde 


twl to-sfpadA to nottsydtl6a 

(7 noreroT ot saad Tetnomniysqxa 

ae 4.) Nottos2-22019 stuart so? etfuesa 
s16Up2 o Noteol Yor sisd (ssnsmiysgxs 

wer as4 . cose eng MOFZIS2-220%9 aysupe Yot 2d iueah = 
164 yYsFuUovrlD yo? edivesh to yvemwe 

ot er ear 59 sis6up2 vot edinesh to viene 


iozirsaqnod rot s2ftwl to SipnA to moryavdtied 
aeeee 8 @ eee ‘nee “ seeee see vee veetese eee eee nexnifosqe 


loziysqnod 10} stad Metnemtyaqx3 


Yebiq2 mov? tefwT To SipnA to nofistusteJ sr.a 
oe, on” avekouh no: Os. are ‘tate wtieeeceee SPNEDGOA OpssToV . 


BO At, de, bad tasT sltensT sot steal Tatnemtveaqx3a, ef.8 


7 Not edfueoR to yvemmue 


oe aa {ues slizast +o seine batemts23 
~~ 7 is 
: A Th a sa 
$id eee weds sf VSluovid to not eo "OF i A 
5 > : a oa Di pae 


all 


+ . - ra 7 
Hk : TILE fA? a2) nc ff: 
tL 21 nie TO ig 


sddebavuesssae FeSt STtenS 






wet eget 


> -F- FSP LFS L/L 


LIST OFs FLGURES 


Page 
SLE ES TCNS LOW pk ee ae oes pag ta ss o's 64 
Typical Rectangular Element in Simple Shear ...... 27 
Some sear (Uy k= sl UU (0) eas trees os 65 
ST UO RAE OMT A 2 gs Coes 66 
Classical Elastic-Plastic Theory Using the 
AUGHMEOCTOSS RALG. (ly Keel ies. cu cc cre eats ee ee es 67 
Hypoelastic Theory using the Jaumann Stress 
ee CCA Re Lion’, 3. cs Meme tints «eis eictans og Hien sian ies 68 
STOnbor Gl sina t Wo abate Unde usbUre BLOVS LON osc cn als 04) 37 
SOUaresCYOSS=SeCLION Of Wie bar cccccnccceeseeceds Sa 
Typical Gridded Region for Numerical Calculation .. 54 
Moment for Torsion of a Square Bar (u/k = 1000) .. 69 
Maximum Stress Function for Torsion of a Square 
Sea (Eyal eee 18 (010 9 ge ie rare ee 70 
Maximum Shear Stress for Torsion of a Square 
pater Ree lU00) oo came. oct er cease ents tes 7| 
Moment for Torsion of a Circular Bar (y/k = 1000). Ve 
Spider Apparatus Mounted on Specimen ............- 74 
Peversntis: andy ralsduccmemrecn scan ccs tra tees st. 75 
Galtbracion- or the Dynamometer’ .....0 ss ccse cece ss 77 
Moment ‘Calibration CUYVemenr ian ares seer cre 102 
Calibration of Angle of Twist for Lower Spider 
Of Ciréular: Specimen) hes he del es 103 


ix 





















fg Leucusacavceecesveee (OO0f = H\u) nofenetxd siqute 


\3 ..e.-. W682 slomte nf tasmel3 yetupnatosh fteorgyt 


~ 
- 
a | 
¢ 
— 
_—) 
f 
= 
— 
i 
~ 
w 
~ 
= 
—_ 
“ 


ait onrell vaosAT sri2zefdesttesia isaieesty 


PO aeeeen etecereewsesce (P= A\y) SPBA Ceaive niet sh ae 


. insiweG siz prfev yyosnt atdeatsoqyt 
5 ead Ue ay omaen woes (CF = S\u) SBR 
\ notetoT sw ysbnU 169 siismeri? Spee 
‘1 sid to sorfas2-22o%9 steupe 
lB. TSU oT 0 rQoH pabbi 40 [62 iqyt 

«y, 


sisupe 6 to nor2%oT ot tasmeM 


. sisupe 6 to nore ot norsonud zesyt2 muinPXeM 
- 32m ee AP ery aor vacees (OOOT = WX) WBa 


c 


| y76upe &. TO NOfexol YOt e2zayvt2 4eade TEM : ¥ 
ry of Lots Cee ness saveiesled nite ieee ee (OOO amma nn 


ui 





&\ teterepeceese MOM D592 No botnvOM eussisqgA yebree . Fa 
hh | v 





» 


s ‘+ , 7 = 
ree ee Te ae evsoubensyT brs amma Yeved Ss 


* oe arr) id oe ii isimennee ai ‘Yo, pores ratte ; ae : 
saoeyeersnyethehaveny Maange naira in 


A i a . 





+ cal 
4 


tend tf 5. 
> > i 


be - Hi 


4 


ve 
4 


Figure 
4.6 


4.7 


4.8 


LIST OF FIGURES (continued) 


Calibration of Angle of Twist for Upper Spider 
OheCERCUloY oPeC Meni mc etre tet cic WN oes a we nes 5 


Calibration of Angle of Twist for Lower Spider 
Areca ae SPECIMEN ici cemimented: Asi slais.d ea eg agentes a 6 


Calibration of Angle of Twist for Upper Spider 
G@psqudre Specimen 4. Beds. Mesae. Siawie.Shoar-..+s 


Moment for Torsion of Circular Specimen .......... 


Experimental Moment for Torsion of Circular 
STSTSGLINEE TS a eee DION reco 2 cots ata anaereea aea arerer rare 


Experimental Moment for Torsion of Square 
SEI Me Nets se cies 4. cca eet ie ere sha t's alen's's Gee ha sss 


Calibration of Angle of Twist for Lower Spider 
GpeCOMmpaL SOMesPECIMGN spire eect 6 ee ose secie we ees 


Calibration of Angle of Twist for Upper Spider 
OP GombaGi SON SDCCIWGM eae cites cer 6c on wate cress oe 5 


Comparison of Angle of Twist Per Unit Length From 


Strain Gauges and Spider Apparatus .........+seee0- 


MHERIENSTLe (e6St SPCCINGN Mere. for. cs sc ce ese es 


Cross-Sectional View of Tensile Specimen SA 
Thesbocationcot. Che sthatneadUCeS .0.... ++ cs sss 


Experimental Tensile Results for a Work-Hardening 
GCA a Pe Semen So 20 Gt Gai eee Seer 


ENG shorsiOl, SPCCIMCNS a wnt eecttete «otk orice ae lereae 
Estimated Slopes of Tensile Results .............. 
Theoretical Curve for Simple Tension Problem ..... 
Simple Shear - Work-Hardening Material ........... 


Torsion of Work-Hardening Material with a Circular 
CYOSS=SOCETON: vee s bie ah tren Sree i oss cn betes aw Oey 


108 


109 


110 


11] 


ie 


85 


86 


lis 

76 
129 
130 
13] 


132 


. 
j 
‘ 7 
(bsunfinos) 23hUdla 40 Tel 
p64 sup ia , 
btq2 tsqgU vot Set] to slonA to nofserdtfsd 4 
Bor PP on re eet namiasqe ysfuott) Fo ana 
‘abfqé 1swol vot tetwT to sinnA to norseydtis2 va 
SPT EAE Se es: o..o. MSmTI9q2 S¥EUpe. To 
Ysbfge ysqqll tot Terwl to SipnA te norsetat {eo BA 
| eee aye. ME .+- HoNnhasqe StBUpe To . 
nontoogé yslusyhD Yo aofedol yo? JrenieM e.# 
i 2 
) isfuovrd to nol2vel yo? InsmnoM -fesasmrisakg Or .B - 
BOI iT Pe eee Pt .«. MoNTOSge a 
aisupe to norero) tot InamoM Tegnetiiveqxg Fra . 
eof ce eo sein a Fae & he tg Se chseis: 8 Bip Ae eel emesis 
vabide 19Wwol.19t J2iwl io slonA to MoNtetdrisy «6h : 
Off oa en eee festuaccds easy ee NSKEGSGe YidZiigaiha) Fo . 
Yabiq2 yaqqU vot Tetwl to Sioa to adtteddifsd ci - 
[rt ere ee be o05%% ome ro at NoMnrssqe “NoeMsgno) to 
nord past tii a9 t2twl 36 sl paA-to Aoetyeqno) BTA 
orl. + sata ... 207616004 Ysbtge Ons zagusa niew2 — 
; £4vashetertee. -...- hentosqe test sftensT-sAT afb 
prtiwede namnragq2 sl rai9i 40 wot). TenéMasg-ecou) 86a 
a8 a isoh 6 Car ee eopuse HrBit2 ent to notis50) sis 
Hnitnebasl-aAvoW 6 yor 2s Ty2o9 al heasT Setnemtysgnd  VF.8 
aD,» Galle wih, ose ei Mats BA bb eee Sat wd eed reese (SEYSTEM 
at Ce ee? eee, roy MTR” , engnitoong woes odt Bf, 





ae sevecveos + ROR afta angele bensiiaas . 


Figure 
5.5 


6.1 


6.2 
643 


LIST OF FIGURES (continued) 


Torsion of Work-Hardening Material with a Square 
CBOSS a OCOU TON .2:/ gi ieeme mene at Wet cpae mkieii Rie’ y's wn a5e os 


Work-Hardening Function for Commercially Pure 
POONA UN ce aaa EE. Gk ad sku a eaid nics Vode os 


Mohr Circle for the Simple Shear Problem ......... 


Typical Element in a Body During Simple Shear .... 


xi 


146 










(baunftnoa) 23valy Se ¥ 


gisub2 6 dtiw fetveisM pntnsbysh- -artoW to sot enol ve 
5 ee oe ee ee sh: Fe “Ou 
sid yifstovemmod ot fot torn” mee “ha 
Me! ss eeakiaeaeer On Petia sicak role Med fun a 


anf cease “eSK2 slami2 patyud Ybos & ni jnsiatd eo tayt 


NOMENCLATURE 


time 
Lame constants 
Young's modulus, Poisson's ratio 


non-negative scalor factors of proportionality 
in the theory of plasticity 


coefficients in the hypoelastic constitutive equation 
yield stress in pure shear 

dimensionless constants, of = aue/ke . Ke = 1/0" 
yield function or plastic potential 

work-hardening function in the theory of plasticity 


work-hardening functions in the theory of hypo- 
elasticity, h = H/2uU 


upper and lower bounds on H 


semi-width of square cross section in torsion problem 
radius of circular cross section in torsion problem 


mesh size in finite difference procedure for square 
cross section 

dimensionless mesh size 

voltage readings taken during experiments 

number of revolutions of the dividing head 


equal intervals in dimensionless deviatoric stress 
during experimental tension test 


estimated error in strains during tension test 


xii 




















TAU TAOWAMOA 
omtt z aa Si 
2inetenod Sm eh 


otiar. 2'doettov , 2ulubod -2  pruory 


| vitisnotdoqerq to 210t26t wohess ay fT apen=non 
vttoties[g 70 yrost ont nf 


_pottsups svrtld tiznoo ottzslaoqyd ont of etnatarttsea 
yoorée Bwuq ar eeavte blaty 

. ¢. . $4 Be | i 
of = "My “ANUS. * To , 2dnsd2znod ezshnofensmts 
fsrtnstoy orteslq yo nortonvT blary 
vi fotdesiq to yiosds sat ni rors: UT pa fnsb46i~ eaee 


“ogy to yoo? sid af enottanut painebtsd= ‘10 
_LS\e = ,vitotsesls 


H no 2ebnued yswol bas ySeqU 


maldoig nohevos af nofsos2 220% svsupe To Atbrwetage 
neldatg nofeios nt noftase 2209 ‘eTuartts to au tDBY 


“i 


Sieupe WOT Siubsoo01g Ssansisttib st tart af asta. team 
nofrsose 22048 


ste degm eeatnotenamtb £; 
bs . 


‘etnemtreqxs pait-vib navet Jepntbsey aps Tov: 
eine 
ad pnkbivib. ais to zhotsufover. %; padi 


vine 












= a 7 







nw ron ey = 


dd 
a 

a 

¥ 


Siraistvet Tp >nomnrb are: 
™ a - . sony +e | P\a 7 
- a ooh ‘ae fino Se on iy: 


= 


truncation error in the slope de/dT 
round-off error in the slope de/dT 


coefficients in the sine series for the slope de/dT 
Vesmlne sss) 


stress tensor components 


deviatoric stress tensor components 


hydrostatic part of the stress tensor 


dimensionless stress tensor components 


dimensionless deviatoric stress tensor components 


velocity vector components 

displacement vector components 

strain tensor components 

strain rate tensor components 

rate of deformation tensor components 

Spin tensor components 

strain increment tensor components 

deviatoric rate of deformation tensor components 


deviatoric strain tensor components 


hydrostatic part of the strain tensor 


hydrostatic part of the rate of deformation tensor 


kronecker delta components 
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NOMENCLATURE (continued) 


In, Il, III, first, second and third invariants respectively of a 
second order tensor Raj 

Be hh ks 

She Chg tomo ink 

ee Ch tom nk 

ee Cet em nn nk invariants used throughout the thesis 

He Deh 

ye 

6 

J = II, 

K = IIT 

W external work per unit volume 

M torque in the torsion problem 

MM, upper and lower bounds on M 

6 angle of twist per unit length in the torsion problem 

8 rate of angular twist per unit length 

E angle of shear in the simple shear problem 

Yow stress functions in the torsion problem, w = ¥/2u 

Y= “2 rate of change of the stress function with respect to 
angle of twist per unit length 

«) warping function in the torsion problem 

S measure of length along the lateral surface in the 
torsion problem 

n external vector normal to a boundary 
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NOMENCLATURE (continued) 


Was components of dimensionless deviatoric stress 

r maximum shear stress 

€ deviatoric strain in the simple tension problem 

S a measure of strain in the simple shear and simple 


tension problems 


Note - For the equations which are in tensor notation, the summation 
convention is implied 
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CHAPTER I 
INTRODUCTION 


The mathematical study of the stresses and strains in plas- 
tically deformed metals began a century ago. In 1864 Tresca proposed 
that a metal yielded plastically when the maximum shear stress attained 
a critical value. A system of five equations governing the stresses 
and strains in two dimensional flow was proposed by Saint-Venant in 
1870 [23]. Saint-Venant postulated that the directions of maximum 
shear strain rate coincided with the directions of maximum shear stress. 
In 1871 Lévy proposed three-dimensional relations between stress and 
rate of plastic strain. In 1913 von Mises suggested a yield criterion 
which is probably the most widely used today. Von Mises also inde- 
pendently proposed equations similar to Levy's. In 1926 Lode conducted 
a series of experiments showing the validity of the Lévy-Mises stress- 
strain relations to a first approximation. In 1930 Reuss, following 
an earlier suggestion by Prandtl, generalized the theory so that elastic 
strains were accounted for. Schmidt (1932) and Odqvist (1933) showed 
in different ways how work-hardening of the metal could be incorporated 
into the Lévy-Mises equations. Little progress was made in the next 
decade. However, since World War II significant advances have been 
made in the solution of special problems. 


A large amount of the difficulty encountered during the solu- 
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tion of elastic-plastic problems arises because there are often two 
different regions separated by the elastic-plastic boundary. Hooke's 
Law is associated with the elastic region while the elastic-plastic 
constitutive equations govern the plastic region. The elastic-plastic 
boundary is found as part of the solution. 

In 1955 Truesdell [1] proposed the theory of hypoelasticity. 


where the constitutive equations are of the form; 
rate of stress = f(stress, rate of deformation), 


where the function f is a homogeneous and linear function with respect 
to the rate of deformation. Chapter II shows the development of the 
hypoelastic theory and shows how the hypoelastic equations may be used 
to represent an elastic-perfectly plastic solid or a work-hardening 
material. The advantage in using hypoelastic theory is that the entire 
body is governed by one constitutive equation and the problem of finding 
an elastic-plastic boundary disappears. 

Truesdell applied hypoelastic theory to Simple problems in- 
volving homogeneous deformation. A.E. Green [14,15] also applied hypo- 
elastic theory to the problems of simple extension, simple shear and 
torsion of a circular bar using a particular form of the hypoelastic 
equations. Solutions for loading were fitted together with solutions 
for unloading. Green obtained solutions using a dynamically correct 
form of the rate of stress. The solutions were then simplified by 


assuming that a certain material constant (u/k) was large with respect 
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2a 
to unity. In 1960 Bernstein [6,7] showed that by changing the initial 
Stresses a particular hypoelastic constitutive equation could describe 
two totally different materials. 

Chapter III contains the solutions of some Simple problems 
using both hypoelastic and elastic-perfectly plastic theory. Compari- 
son of the two theories is made wherever possible. It is shown in 
Chapter III that under certain conditions the dynamically correct form 
of stress rate may be replaced by the material rate of stress and the 
solutions obtained with this approximation are the same as the simplified 
solutions obtained by Green [14,15] in a different way. This is important 
because the use of the material stress rate simplifies most problems 
considerably. This thesis applied hypoelastic theory to a non-homogeneous 
deformation problem (torsion of a Square bar) and uses the material 
rate of stress. 

Chapter IV describes some experimental work that was carried 
out using a work-hardening material. Torsion and tension tests were 
performed on commercially pure Aluminum specimens. 

In Chapter V various work-hardening problems are solved 
using a particular form of the hypoelastic constitutive equation. The 
work-hardening function is evaluated from the results of the experi- 
mental tension test. This function is then used in the hypoelastic 
constitutive equation to predict the results of the torsion experiments 


for a square and a circular cross-section. 
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CHAPTER II 
THE RELATION BETWEEN HYPOELASTICITY AND PLASTICITY 


2.1 The General Hypoelastic Constitutive Equation 


The classical linear theory of elasticity is based on the 


relation 


stress = g(small strain from an unstressed state) (2.1) 


where g is a single valued function. This relation has been generalized 


for the theory of finite strain as 

stress = h(finite strain) (2.2) 
where h is a single valued function. Truesdell [1] suggested that the 
basic assumption (2.1) of the classical linear theory could be regarded 
as a first approximation to a more general theory of finite strain dif- 
ferent from (2.2) and proposed the relation 


stress increment = f, (strain increment) 


or 


rate of stress = Ff, (rate of deformation) 
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where f, and f, are single valued functions. 

Further, Truesdell sought to free this concept of elasticity 
from any connection with a natural or otherwise preferred state. It 
was proposed that the constitutive equation may depend on the current 


stress, so that 
rate of stress = f(stress, rate of deformation) . (2:33) 


This theory may be regarded as representing a material in 
which the stress is built up by increments which at any instant obey 
a linear stress increment-strain increment relation with coefficients 
depending upon the current stress. The response predicted by equation 
(2.3) is identical to that predicted by the classical theory if the 
strain is infinitesimal. 


Equation (2.3) can be written as 


k 
ee a Ue (2.4) 
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where 


) 


Sere 
13J J31 


is the rate of deformation tensor, Vv, are the velocity components and 


the semi-colon denotes covariant differentiation. The stress tensor 
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is denoted by oF while Doy/Dt is the "rate of stress" which may take 
one of various forms each of which must be objective in order that the 
constitutive equation be dynamically correct. 

The rate of stress that is used in this thesis is the co- 


rotational derivative, 
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are the components of the spin tensor and 
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is the material derivative of the stress tensor. 5 

Prager [2] has discussed certain other definitions of stress 
rate and has noted that the definition given by equation (2.5) has the 
property that zero stress rate implies that the stress invariants are 
stationary whereas the other definitions do not have this property. 

The stress rate (2.5) used in this thesis was derived by 
Noll [3] and earlier by Jaumann [4] in a different way. The Jaumann 
stress rate measures the material rate of change of the stress components 


with respect to coordinate axes that rotate with the material and in- 


a 











Pe. eee . ‘ tite negl Galas 


ath Ls ‘y th dood état Sere te 


brain t ta ‘iy J ene 


e"Su7L ae svn vP9at lufvetaa nd 2 
~~ 


) Korcageth ca ee] age 

quvtp wel ATS Gt Sek rir sar oo 

ery) naiwavnped add: Fe43 i fvge i anti al2 ores Jan} i 

Maegan ade wuesditen ‘iti adet Tats ionsit SH rainy Vv PaaGt 
i, i'n bowl aM 2 Prams fi fi teoy. {2.2)-$009 24913 2-907. 


4 wat mie a i 
fil ini 6° nn) ie ; | 
e + a Pe ae + - 7 - mo bia fn 


s 


av 








Fim . 
- - 
- a 


17 





Stantaneously coincide with the fixed axes. 
Returning now to the general form of the constitutive equation 


for hypoelasticity, equation (2.4) can be rewritten as 


where 


are the components of the stress tensor with twice the elastic shear 
modulus as the unit of stress. 

All constitutive equations must satisfy certain invariance 
requirements. In order that the well known principle of objectivity 
be satisfied, it can be shown that the function f must be a hemitropic 


function of d and T, where f, d and T are the matrix equivalents of 


ek m : 
Fas d, and t, respectively. 


According to a theorem of Rivlin [5], a hemitropic polynomial 


f of two tensor variables d and T can be represented in the form 


Tuan (2.7) 


2 
I + a.,d + a,,d° + Any! 02! 


F(d,T) = agol + ajod + agod 


Td + dT) + ay 9(aT? i Td) 





*For simplicity matrix notation, which is valid in general coordinates, 
is now used. 
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Use has been made here of the Cayley-Hamilton theorem which states that 
a square matrix satisfies its own characteristic equation. That is 
a 2 

= 1,7 - I1,T - IIT, 
where I, II and III are the first, second and third invariants re- 
spectively. The coefficients Avg above are functions of the ten in- 


variants Ty: Il), IIT, I, Il, IIl,, heave sa0-OT which the last 


four are the joint invariants defined as 


ty Jka _ ik 24m 
», Wise e BARS 


By using the principle of material objectivity, the function 
f has been put in the form of equation (2.7). This equation can be re- 
fined further by applying dimensional invariance. 

One of the hypoelastic hypotheses is that no constitutive co- 
efficient of a hypoelastic material can have a dimension independent of 
Stress. Since there is no modulus with dimension "time", there is no 
relaxation effect. 

Since pti /ot and ds are of dimension tl where t is time, 


* 
equation (2.6) must reduce to the form 


ee ee ee ee ee 
*The more general bodies which do not have the restriction of linearity 
in dt are called hygrosteric by Noll [3]. 
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pe Ay, Sy (2.8) 
where A is a dimensionless function of tie It may be deduced from 
equation (2.7) that og = Ap] = Ano = 0, Ang? 8o) and Ano are of degree 
me ik ‘ k hae 2 
one in do and Aig? 97] and aj5 are independent of qd, - The joint in- 


variants 


Kk 24myn 


m 3 
d, 7 and a0. = t, monk 
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Since ag are power series in T, equations (2.6) and (2.7) 


do not appear since they are of degree two in d 


reduce to 


ara § Oty I + ad + 6 oT + Rat (2.9) 
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for a compressible solid, where 0, (y = 0,1,2,...,11) are dimensionless 


analytic functions of the three principle invariants of T: 
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= A pu + (a, + F)d + ag lat + Td) + 
(aR + aN) I + (aR + Oy NT + 
(aig + ay,N)T> + Fag(at® + Td) 
If 


a, (00,0) = 1 and at) (00,0) RS 


(2.10) 


where A and pu are Lamé constants and v is Poisson's ratio, the constitu- 


tive equations of hypoelasticity reduce to those of the classical theory 


of elasticity as a first approximation. 


However, the constitutive differential equations (2.9, 2.10) 


and the initial conditions are not enough to define uniquely the hypo- 


elastic material behavior. A hypoelastic material is defined by the 


assignment of a set of constitutive differential equations (2.8) and a 
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corresponding equivalence class of stress configurations. Bernstein [6,7] 
has shown that by changing the initial stress, a particular constitu- 


tive equation may define totally different materials. 


2.2 An Elastic-Perfectly Plastic Constitutive Equation 


The following discussion assumes that the elastic part of the 
strain is small and can be regarded as infinitesimal. The rate of de- 
formation tensor qi may then be written as the sum of an elastic 


part di and a plastic part di [31] 


d; 32 des + dP. 
ig 1J iJ 

At this stage it is convenient to consider strain increments rather 

than rates of deformation so that the total strain increment de. « can 


J 
be written as the sum of an elastic and a plastic part. 


de.. = de® + de? 


ij ‘4 ij (201) 


The strain increment de. [9] is defined by the equation 


where du. is the incremental displacement of a point whose current 
position-vector is x referred to fixed Cartesian axes {Ox}. It 


should be noted that du, does not denote us (x + dx) - u. (x). 
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If the total strain is small the total strain rate is identical 


to the rate of deformation and can be written as 


where Vv. are the velocity components referred to rectangular Cartesian 
axes {Ox, } and a comma preceding a subscript denotes partial differen- 
tiation with respect to the appropriate spatial variable. 

An assumption of classical plasticity is that there exists 


a yield condition described by the equation 
Flo; 5) = 0 


where C5; is the symmetric stress tensor and the sign of the yield 
function f is chosen so that f < 0 corresponds to the elastic regime. 
The function f also depends on the temperature and prior history of 
work-hardening and heat treatment. 

The flow rule [8] for perfectly plastic materials associated 


with the yield function Flo, 5) is given by 


Reno Peet 
dis r ip (2512) 





where it is assumed that the yield function f is regular at the stress 


point considered and di are the plastic rates of deformation. i is a 
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non-negative scalar factor of proportionality and is not a material 


constant, but varies during the deformation and 


MeeaCerit +f <0) -oremetmeia Ole and << 





h>O if f=0 andf=0. 


The function f is often called the plastic potential. An example of 
such a function is the von Mises yield function, which is considered 
later. 

For a perfectly plastic solid with a constant state of stress 
at the yield point non-vanishing strain rates may occur. However, for 
a work-hardening material, an increase of stress beyond the yield limit 
is required to produce non-vanishing strain rates. The flow rule for 
a work-hardening material with a regular yield surface is 


Ossi fugiest0d one fol toni Qgandeti=n 
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where the positive work-hardening function g may depend on the stress 
and strain histories, but is independent of stress rates. Therefore 
the plastic rates of deformation for loading (f > 0) are assumed to be 


linear functions of the stress rates. 
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It is assumed that metals are plastically incompressible. 


Consequently 


ae 

dik 0 e 
With the assumptions of isotropy and absence of Bauschinger effect, 
along with the experimental fact that yielding is to a first approxi- 
mation unaffected by moderate hydrostatic tension or pressure, it can 


be shown [9] that the yield function can be put in the form 


I 
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where 
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and f is an even function of Il. 44 is the deviatoric stress tensor. 


Von Mises suggested the following yield function for a non- 


work-hardening material; 


where k is a constant. That is 
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Equation (2.12) then yields 


Dede 
diy r 344 : 


In terms of strain increments this can be written as 


Pegs 
des; di Sij ; 
From Hooke's Law 
Hee ote oi do 


ij ee 


where do = doy ,/3 . Equation (2.11) then implies that 
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(2513) 


These equations, which are known as the Prandtl-Reuss equations, may 


be rewritten as 


(a) 
dS, 
tie ee 
(b) 
7g Ce: 0 
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(c) 
Si oS een Le 


(2.14) 
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where 


Equation (2.14b) is the compressibility condition. The external work 


dW per unit volume done on the element during the strain increment 


dW = dW, + dW 
where 


= e = p 2.15 
dW, = 9;5 dey; » dW, =S., de. . (2,15) 
dW, is the elastic recoverable energy while dW, is the plastic work per 
unit volume. Now since plastic distortion is irreversible in the thermo- 


dynamic sense, 


dW, 20 
* 
during plastic deformation. From equations (2.13) and (2.15) it follows 


that 


ee ee ee ee eee, 
*Yield surfaces in stress space which do not enclose the origin are 
not considered. 
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dW, = akm"dy > 0. 


Now during loading in the plastic region 
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Therefore during loading in the plastic regime 
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For unloading or during loading when Saj Sig < ake 
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The above is one example of the constitutive equations used in the 
solution of elastic-perfectly plastic problems. It should be noted 
that the stress increments in the above equations are referred to 
fixed axes and are not objective. However, if the rotations are 


small these equations are valid. 


2.3 A Hypoelastic Constitutive Equation Similar to the Prandtl-Reuss 
Equation 


In this section, it will be shown that under the proper as- 
sumptions the constitutive equation for a hypoelastic body can be re- 
duced to a form resembling that of the elastic-perfectly plastic body 
described in Section 2.2. 

From Section 2.1 it is clear that in developing the theory 
of hypoelasticity, no assumptions are made as to the magnitudes of the 
deformations. The theory is completely general and applies to large 
as well as small deformations. If the dynamically correct form of the 
Prandtl-Reuss equations is used, the only assumption concerning magni- 
tudes of deformations is that the elastic strain is infinitesimal. 

The constitutive equation (2.16) for the elastic-perfectly 
plastic body obeying the von Mises yield criterion is a combination of 


classical linear elasticity and pure plastic flow in which the components 
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of deviatoric strain rate are proportional to the components of de- 
viatoric stress. 

The constitutive equations (2.9) and (2.10) for the com- 
pressible and incompressible hypoelastic body are general in that 
they may describe many different materials depending on what values 
are given for the coefficients in the equations. 

Hypoelastic and plastic constitutive equations are rate 
independent. Consequently the hypoelastic theory can not, for example, 
describe a viscous fluid. Also in both the hypoelastic theory and 
the plastic theory as described in Section 2.2 it is assumed that the 
materials are isotropic. 

Thomas [10] has obtained the dynamically correct form of the 
Prandtl-Reuss equations in terms of stress rates rather than stress 
increments. These equations, referred to a fixed system of curvilinear 


coordinates {oy'}, are 
= Ce ae (aig) 
The symbols used in the equation are defined above. 


Now for loading in the plastic domain 
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Substitution into equation (2.17) yields 
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where the second term on the right is present only for loading in the 


plastic zone. For an elastically compressible material 


qe 
1 Do _ Sa ckk 
aeitaas (1 + eS (2.19) 


Consider, now, the general constitutive equation (2.9) for 


a compressible hypoelastic material. 
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then equation (2.20) becomes 
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Equation (2.22) can be written in terms of the deviatoric and hydro- 
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etgt, = Ske ke (2203) 


Boer (tga. (2.24) 


Equation (2.23) is similar to the Prandtl-Reuss equation (2.18) ex- 
cept that the second term on the right hand side of equation (2.18) 
is present only for loading in the plastic zone while the corresponding 
term in equation (2.23) is present for all "loading" deformations. 
Also in the elastic-perfectly plastic theory, di; is divided into elastic 
and plastic components whereas this is not true in hypoelasticity. Thus 
the second term on the right hand side of (2.18) is not quite the same 
as the corresponding term in equation (2.23) even thoygh we may be 
loading in the plastic region. 

Upon comparing equations (2.19) and (2.24) it is noticed 
that in equation (2.24) "d" is the total rate of dilatation while in 
the elastic-perfectly plastic theory it is assumed that the body is 
plastically incompressible. 

Truesdell [12] has shown that for certain incompressible 


hypoelastic solids, the constitutive equation can be written as 
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The differences between this equation and the corresponding one in 
classical plasticity are identical to those stated above for the com- 
pressible body. 

Further comparison between the two theories may be made in 
regard to the notion of a yield criterion. In the theory of plasticity 
a yield criterionis assumed. Truesdell [13] on the other hand has 
shown that by using the equations of hypoelasticity one can obtain 
what is called a hypoelastic yield as one of the results in the solu- 
tion of the simple shear problem. In the paper referred to above, 
Truesdell solves the problem of simple shear using equation (2.23) 
and compares the hypoelastic yield to the von Mises yield for various 


values of (u/k). 


2.4 A Hypoelastic Constitutive Equation Suitable for Some Work- 
Hardening Materials 


The hypoelastic constitutive equation derived earlier in 


this chapter can be written as 


Dees 
Se 
Ot de +h Tee deel ys (2.25) 
where Lr = $4 4/4u and h = - aue/ke. However, in general h is a di- 


mensionless function of J and K where 
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In Chapter V, this thesis considers the application of 
equation (2.25) to problems involving work-hardening solids. The 
assumption is made that h is a function of J only. The functional 
relationship for commercially pure aluminum is approximated from the 
results of a simple tension experiment. The constitutive equation 
(2.25) is then used to solve the torsion problem for this material 
for both square and circular cross-sections, and the solutions are 


compared to experimental results. 
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CHAPTER III 
SOLUTIONS TO SOME HYPOELASTIC PROBLEMS AND COMPARISON 
TO CORRESPONDING ELASTIC-PERFECTLY PLASTIC SOLUTIONS 


This chapter contains the solutions of some problems using 


the hypoelastic constitutive equation 


We d!. - of M! 
pe ig ON Tay 
Lie ee ti . 1 : 
where a = 2u'/k” and M' = They og: Comparisons are made to the 


classical elastic-plastic solutions in which the Prandtl-Reuss equations 


and the von Mises yield criterion are used. 


3.1 Simple Extension 


Consider the velocity field 
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for an incompressible material where the components of velocity are 


referred to a Cartesian coordinate system. Then 
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and the rate of spin tensor is zero. Assume 
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The constitutive equation can be written as 
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Now since Wij = 0 and the stresses are assumed to be homogeneous, 


equation (3.1) becomes 


ela 


(372) 


Equation (3.2) and the initial condition that T = 0 when c = O implies 


[14] that 


T = +e tanh (= ac) 


25 






attinw ad neo notisups sv roses (senoa or 


' 
1d , 
= - if 
; > 1} i! -. j | i 
Lie ; ‘eyMoy! . eae tar ar sU 


, 2vosnepc nod od ot bomuzes ave eseesvte Sad bos 0 = som til 


uf 
(f.€) 8 natden 


2smoosd 


4ayomevsg onfesoiont yi fsatnotonom 6 8d od bemwers 


26 


or 


1S E> 


) tanh (Y3Z=c) . 


A{e 


Ta 
T= > | 


The classical elastic-plastic theory predicts that 


and 


where "c" is the logarithmic strain. 
Figure 3.1 shows the relationship graphically between di- 
mensionless deviatoric stress and logarithmic strain for both the 


elastic-plastic theory and the hypoelastic theory when u/k = 1000. 


3.2 Simple Shear 


In this section the problem of simple shear is considered 
for both the elastic-plastic solid and the hypoelastic solid. The 


problem is solved for both cases using the objective stress rate + 


and the non-objective stress rate : 


The velocity field, referred to fixed Cartesian axes Ox. is 


Vie 2c(t) Xo » Vo = V3 = 0 
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n ac 
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CO 7 
th 


and 2c = tan € where € is the angle of shear as shown in Figure 3.2. 


The equations of motion are satisfied if body forces are 


zero and c = Q. 





Fig. 3.2 Typical Rectangular Element in Simple Shear 
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The rate of deformation tensor is 
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while the spin tensor is given by 
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where 343 is the deviatoric stress tensor and 
T= Tt e=s5(t) 


That is, it is assumed that a state of homogeneous stress exists. The 


deviatoric stress is used rather than the non-dimensional deviatoric 
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stress in the simple shear problem. 

The elastic-perfectly plastic solutions are now obtained 
using both the dynamically correct form of the Prandtl-Reuss equations 
(i.e. with the Jaumann stress rate) and the approximate form with the 
material derivative of stress. The hypoelastic shear problem is then 
solved using both the Jaumann and the material stress rate. 

If the Jaumann stress rate is used, the elastic-plastic 


theory predicts that 


See 
peed J = 
Dt * Sak? ky + SjK ei = 2H 955 
during elastic deformation. It follows that 
dig, = 
Ar cs (se) 
and 
ds 
dou =a ce eit : (3.4) 


Integration of equations (3.3) and (3.4) subject to the initial condi- 
tions that T = 0 and S = O when c = O leads to 


T = u(1 - cos2c) 


29 


n 


, 
hs a 7 













oe : i _ 
7 Leite 

neidorg "sare shqmte od af ezewe 

: 7 

¢ anortuloe aftesiq yitosisq-ottesis off ; 


snattsuoe e2uan-Tst qd sdt to mrot tosii0o VITSSimaNND Sas tod 1 pate 
‘or otemrxorags 9dt bas (etsy eesyd2e ansmuet + ais nittw st r) 


~ j | > F > « 
um ob matdaa sti2 otteslseqyd eit <e2arte Yo evitey inab | isivssem 
: ol > rye ‘ i 
stystsm ant bas nnsmusl snd rfod pateu /10¢@ 


oe 


fosoit2el + boou ef Stat 2eeavwdse finest say Ti 


Js 239 rbsy ¥ 





Tae otal ole a ta 
i fertint ant ot t9etdue (B.€) bis (E. €) anotssips to norter 
bs 


ov as e 
od abaot 0 = 9 ast O.='2 bns 0 = 7 ddl 


hie oe heal Sadie 









t ; e 





and 


See uesiteec. . 


The yield criteria predicts that the material yields when 


2 
cos2c = 1 - _ 
2u 
Then 
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For the yield point to exist equation (3.5) implies that 
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The Prandtl-Reuss equation 
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during plastic deformation. In other words the maximum shear stress 
remains constant. This is merely a statement of the yield criterion. 
To find the individual components of stress during plastic 


flow, the equations 


1/2 
dT i AAG: 
T= 2(1 - 31) (k2-T7) 
and 
1/2 
S = (k* - T°) 
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are solved subject to the condition 
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If the material stress rate is used in place of the Jaumann 


derivative then 


T= 0 
éC as 
= [ui 
S’ = eC 
and 
T=0 
Ze > kK 
Say Ul 
S =k 


For most metals k/u << 1. Consequently these two solutions 
are almost identical. 


The hypoelastic constitutive equation 
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in which the Jaumann stress rate is used, reduces to 
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for simple shear. Again the initial conditions are 


(S) <9 = 0 
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For S to remain real it can be shown that 
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If the material rate of stress is used, equation (3.8) reduces to 
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and 


subject to the initial conditions 


(S) <0 mt 
and 
(T) 20 “ 
Therefore 
T=0 
and 


Ss) ettanh (As) 


The results of the four solutions are shown graphically in 


Figures 3.3 to 3.6. Figure 3.3 shows a plot of t/k vs. tan & for various 
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values of u/k where 


and 
Cantera 2C-, 


For the values of u/k used in Figure 3.3 the results obtained using 
the material derivative are the same as those obtained when using the 
Jaumann derivative. It can be seen that in the hypoelastic solution, 
t/k rapidly approaches unity from below. 

Figure 3.4 indicates however that for u/k = 1 the same re- 
sults are not obtained when the material stress rate is used in place 
of the Jaumann stress rate. In fact a a erevent yield point is ob- 
tained when the material derivative replaces the Jaumann derivative 
in the elastic-plastic theory. Figures 3.5 and 3.6 show the values 
of T/k and S/k for u/k = 1 when the Jaumann stress fa tests used. 

It can be shown that if the Jaumann stress rate is used, 
T/k approaches k/u asymptotically from below for both the elastic- 
plastic and the hypoelastic solutions. At the same time S/k rises 
above ¥ ] - Ke pu? and then approaches v1 - Kent asymptotically from 
above. The only difference between the elastic-plastic and the hypo- 


elastic solutions is that t/k = 1 in the plastic solution whereas t/k 
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approaches 1 from below in the hypoelastic solution. 

Two conclusions can be drawn from the above exercise. It 
appears reasonable to conclude that hypoelastic theory can be used as 
an approximation for elastic-plastic theory in certain problems. This 
is beneficial since in some cases the hypoelastic constitutive equations 
are easier to solve. Also since p/k is approximately 10° for most 
metals, the Jaumann stress rate can be replaced by the material stress 
rate if the components of the rate of spin tensor are of the same order 


of magnitude as those of the rate of deformation tensor. 


3.3 Torsion 
3.3.a Statement of the Problem 

Consider a prismatic bar of any cross section twisted by 
couples applied at the ends. The lateral surface of the bar is to be 
free from external forces while the ends are allowed to warp in the 
axial direction. 

The X3 axis is chosen parallel to the generators with the 
origin at one end as shown in Figure 3.7. 
3.3.b Classical Elastic Torsion 

Saint-Venant [17] solved the classical elastic torsion problem 
by the semi-inverse method. That is certain assumptions were made con- 
cerning the deformations of the twisted bar and it was shown that a 
solution could then be found that would satisfy the governing equations 
and the boundary conditions. Then from the uniqueness of elastic solu- 


tions it is known that this is the exact solution. 
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Fig. 3.7 Right Prismatic Bar Under Pure Torsion 


The displacements are assumed to be 


where 8 is the angle of twist per unit length and $(X] »X_) is called 
the warping function. Displacements and rotations are assumed to be 


small so that the boundary conditions may be referred to the undeformed 


configuration. 
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Now since the strain tensor is defined as 


it follows that 
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From Hooke's law it follows that 
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where py is the modulus of rigidity. If the stresses are to satisfy the 
equilibrium equations and the boundary conditions then it can be shown 


that 
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on the boundary, where S is a measure of length along the lateral surface 
in the X1 Xo plane. Body forces are neglected in the equilibrium equations. 
The problem can also be formulated in terms of the Prandt] 


stress function [18]. If 


ee 
Cia aX (X7 2X5) 
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where ¥ is the stress function, the governing equation is 





Subject to 
Y= constant 


on the boundary. For a simply connected region ¥ will be chosen to 
be zero along the boundary. 


It can also be shown that 
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where M represents the torque applied at the ends of the bar and the 
integral is taken over the area of the cross section. 
3.3.c Hypoelastic Torsion 

Consider a hypoelastic material that obeys a constitutive 


equation of the type 


i pe di 5 - Seo Si (3.9) 
This equation was developed in Section 2.3 for both compressible and 
incompressible bodies. 

Throughout this discussion Cartesian coordinates are used 
and a suffix preceded by a comma denotes partial differentiation with 
respect to the appropriate spatial variable. 

Displacements and rotations are assumed small so that the 
boundary conditions may be referred to the initial configuration. 


The velocity field is assumed to be 


V3 = 8 (x) 9X5 96) 


where 6 is the rate of change of angle of twist per unit length. Then 
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Upon replacing the Jaumann stress rate by the material stress 


rate, equation (3.9) becomes 
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last two terms in the above expression are omitted and the material 


derivative can be replaced by the partial derivative with respect to 


time. The constitutive equation (3.10) can then be written as 


oT. 
a a ee Labor 


gh 


f J mi? 7 
f |B 
aa — - ees 
xed $6 +0 
ot 2eoubsy (17.8) wervee 
es 
16 sn 10 
! a 4 + — ~™ 4 pL = ats aon 
| pry t a6 3G 







Si? .x to Jnabnsasbnt ad oF bomvzes Sis eater ge sents bts Bis. sont 
i. i +! =e a 
fefyatem siz bos bettima 915 notzesiqxs syade sit nf emied Gar # 


2) 





as toaqee idtw sviteyitab fstiosq and vd bexslqar sam 


26 asttiww od ners nso (Of,£) gots s | rnd Fei 02 3 T tye: 
Ap . - : aa _ a 


_ 


43 


Since 6(t) is a monotically increasing function of time one can write 
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From the initial conditions and the first four of equations (3.12) 


it can be shown that 


for all 6. 


If body forces and acceleration effects are neglected, the 


only non-trivial equation of motion reduces to 








The boundary conditions require that w be constant on the boundary and 
Since the section is simply connected the constant can be taken as zero. 


The hypoelastic torsion problem reduces then to solving the equations 
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where 
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on the boundary of the cross-section. Multiply equation (3.14a) 


by abe and equation (3.14b) by (5 - aT.) and add the two equations. 


The result is 
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(b) 
7- A TT3 #0 
(3.17) 
(c) 
(d) 
z- 8 (Thy + Th) #0 


Multiplying equation (3.12e) by 13 and equation (3.12f) by T53 and 


adding gives 


Pag (M13 + Tog) = Wily - a (TS, + T5)1 
Since W' is always positive it follows that 
7° a°(T%, + Tee le> 0 (3.18) 
and a” (TS, + i>) approaches 5 asymptotically from below. It is in- 


teresting to note that the von Mises yield criteria is 
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7° a°(T%, + 198) > 0 
then 
Z- A173 > 0 
and 
77 ale eat) 


Showing the validity of assumptions (3.17b) and (3.17c). It is necessary 
that 113 # 0 and T53 # 0 for assumption (3.17a) to be true. The as- 


Sumption that T)3 # 0 and T53 # O was also used in obtaining equation 
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Equations (3.13) and (3.21) then become 
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It can be shown that 


Therefore equation (3.22) is elliptic until infinite deformation at 
which time it becomes parabolic. 
The boundary condition associated with equation (3.22) be- 


comes 


oy _ 
ae 0 (2223) 


on the boundary. 

The hypoelastic torsion problem could also have been formu- 
lated in terms of the warping function 9, however it was found that 
problems arose when trying to satisfy the boundary conditions in the 
numerical procedure. 

3.3.d Numerical Procedure for a Square Cross-Section 

The governing equation (3.22) for hypoelastic torsion and 
the associated boundary condition (3.23) will now be applied to a 
square cross-section. Since the classical elastic and plastic solu- 
tions possess eight fold symmetry, it is reasonable to assume that the 
hypoelastic solution will behave similarly. Referring to Figure 3.8, 


only the shaded area OMN in the first quadrant will be considered. 
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Fig. 3.8 Square Cross-Section of the Bar 


In applying the boundary condition (3.23) to the triangle 
OMN it is clear that 
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where n is the normal to the boundary ON. 


Equation (3.22) becomes 


2 

oY ony oY oY oY 
A + 2B +C +D—+t+E=—+F=0 (3.24) 
ax,” 9X4 9Xo ax5° OX) OX5 


Since 1,3 and T53 are zero when 0 = 0, the coefficients A, B, C, D, E 
and F can be evaluated at any point in the region OMN. Equation (3.24) 
is solved for oy/36 using a finite difference procedure. The stress 
function y is known to be zero initially and by using a forward inte- 


gration of the type 
p(6+A6) = w(6-As) + 2A6 aye) 0((a8)%) F 
or neglecting higher order terms 


v(e+Ae) = y(6-Ae) + 2A6 = ayo) : 
~(6+A8) can be found at each point of the finite difference mesh. New 
stresses can be found from the new stress function w(6+A6@) which are 


then used to calculate new coefficients A, B, C, D, E and F. This pro- 
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cedure is repeated until the total angle of twist per unit length 
reaches the desired amount. 
An 11 x 11 finite difference mesh was set up in the octant 


OMN as in Figure 3.9. 





Fig. 3.9 Typical Gridded Region for Numerical Calculation 
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The truncation error in each of the above approximations is 0(a7A*). 
Substituting the above expressions for the derivatives of Y into equation 
(3.24) results in fifty-five linear equations for fifty-five unknowns. The 
fifty-five unknowns are the values of Y for the fifty-five points in 


the octant OMN excluding the eleven points along the boundary MN where 


Y is known to be zero from the boundary conditions. The equations are 
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Each time that the above equations were solved for the dimensionless 
quantities Le (fa new stresses and the corresponding moment were cal- 


culated by the following equations. Since 
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From equations (3.13) it is known that 


These stresses can be found approximately as 
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dx dx. 
dX, = — dX, = — 


] a 2 a 

The double integral in equation (3.25) was evaluated by Simpson's 
Rule. 

The above numerical procedure was carried out using a model 
360/67 IBM computer. The results are tabulated and the corresponding 
graphs are given in Section 3.3.e. 
3.3.e Results for Torsion of Square Bar 

The classical elastic theory for torsion of a square bar [19] 


predicts that 
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and the bar begins to yield when 
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It can also be shown that the perfectly-plastic moment My is given by 
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so that 


and the classical elastic-plastic theory predicts initial plastic flow 


when 


ware 
(6a) = 7357 

The numerical procedure described in Section 3.3.d calculates M/2ua®, 
¥/2ua ,* S43/2u, S53/2u and t/2u** for values of (6a) up to 3(8a) Some 
of these results are shown in Table 3.1 for an increment size of 

ee) os abl 
A(@a) = T50 (7355) . The computer printed out results for the first 
increment, which is merely the classical elastic problem, and for every 
twenty-fifth increment thereafter. 

Graphs of M/2ua> Vs 6a, ¥__./2ua Vs 6a and T1a/2h Vs 6a for 


max 
] 


u/k = 1000 and (6a) = (gem) (qgu5) are shown in Figures 3.10, 3.11, 


and 3.12 respectively. 


3.3. Hypoelastic Torsion of a Circular Bar 


For the case of a bar with a circular cross-section it can be 
shown [14] that integration of equation (3.9) yields 
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pe = x (S) tanh (E re) (3.26) 


where 6 is the angle of twist per unit length and r is a measure of 
radial distance. The Jaumann derivative in equation (3.9) was replaced 
by the material derivative in obtaining equation (3.26). 


The torque is 





a 
M = 2nk | ré tanh (7 re) dr 
0 
or in dimensionless form 
ee? 

Me Ky ee Ushee) ai tt 
mea MG) | OG) tanh CB oadatg) 

) 


where "a" is the radius of the circular bar. 

Simpson's Rule is used to evaluate the torque for various 
values of 6a. For the computations, u/k = 1000 and n = 20. The 
results are plotted in Figure 3.13 and compared with the elastic- 
plastic and perfectly plastic solutions. 


For the elastic-plastic solution, 
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TABLE 3.1 


Results for Elastic-Plastic Torsion of a Square Bar 








(ea) A a Tmax 
Qua pa 2u 
(x10*9) (x10*3) (x10*3) (x10?) 
0.0000 0.0000 0.0000 0.0000 
0.001644 0.001844 0.0004834 0.001028 
0.1233 0.1378 0.03614 0.07663 
0.2465 0.2722 0.07162 0.1504 
0.3698 0.4004 0.1058 0.2188 
0.4931 0.5198 0.1383 0.2798 
0.6164 0.6288 0.1686 0.3324 
0.7396 0.7264 0.1965 0.3759 
0.8629 0.8123 0.2220 0.4108 
0.9862 0.8870 0.2449 0.4376 
1.1095 0.9514 0.2655 0.4577 
1.292/ 1.0056 0.2839 0.4721 
1.3560 130516 0.3003 0.4821 
1.4793 1.0903 0.3148 0.4888 
1.6026 Lel228 Q73276 0.4932 
1.7258 Ts501 0.339] 0.4960 
1.8491 hl / 34 0.3493 0.4978 
1.9724 1.1924 0.3584 0.4989 
2.0957 1.2088 0.3665 0.4996 
2.2189 Lecce, 0.3739 0.5000 
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CHAPTER IV 
EXPERIMENTAL CONSIDERATION OF A WORK-HARDENING MATERIAL 


Two separate tests, torsion and simple extension, are carried 
out on commercially-pure Aluminum. The purpose of the torsion experi- 
ments is to obtain a curve of moment versus angle of twist per unit 
length for a work-hardening material. The tension test is performed 
in order to obtain the work-hardening function which is used in calcu- 
lating a moment versus angle of twist per unit length curve theoretically. 


The theoretical and experimental results for torsion are then compared. 


4.1 Torsion Tests 


4.1.1 Experimental Apparatus 


Since it was necessary to measure angles of twist as small 


as 107° 


rad. an apparatus which would measure small angles of twist was 
designed. The device measures tangential displacements by a linear 
transducer and the displacements are converted to angles of twist by a 
calibration curve. 

As is shown in Figure 4.1, circular mounts are attached to 
the specimen at each end of the gauge length by means of screws. Four 
indentations are made at each end of the gauge length at 90° from each 
other and the points of the screws fit into these. This is to ensure 


that the planes of the "spiders" (circular mounts) are perpendicular 


to the axis of the specimen and also that the outer circular edges of 
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Fig. 4.1 Spider Apparatus Mounted on Specimen 


the spiders are concentric about the axis of the specimen. To the outer 
edge of each spider is attached a very fine wire which is used to 
measure tangential displacement. The displacement is amplified by a 
lever mechanism with a mechanical advantage of about ten before being 
attached to the core of the transducer. Figure 4.2 shows the lever 
arms and the transducers. A displacement of the transducer cores causes 
a change in voltage which is measured by voltmeters connected to the 
transducers. 

Torques are measured by means of a set of four strain gauges 
mounted at 45° on the shaft leading to the upper set of jaws on the 
torsion machine. The dynamometer is explained in more detail in Section 
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Fig. 4.2 Lever Arms and Transducers 


4.1.2 The Test Specimens 


The material used was commercially-pure Aluminum which was 
annealed by the following procedure. The machined specimens were heated 
gradually in an oven from room temperature to 700°F over a period of 
forty-five minutes. The temperature was controlled at 700°F for thirty 
minutes and the specimens were then allowed to cool to room temperature 
overnight. 

Two different cross-sections were used in the torsion tests. 
The specifications are shown in Figure 4.18. The taper was kept to 


within + 0.001". 
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Fig. 4.18 The Torsion Specimens 
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4.1.3 Calibration and Experimental Procedure 
4.1.3.a Calibration of the Dynamometer 


The dynamometer was calibrated by mounting a lever arm on the 
upper jaws of the torsion apparatus and applying known moments by means 
of pulleys and various weights as shown in Figure 4.3. As was mentioned 
in Section 4.1.1, strain gauges were mounted on the shaft leading to 
the upper jaws and the moments were correlated to the strain gauge 
readings which were taken from a voltmeter. 

Care was taken to minimize bending or axial loading during 


the calibration. However, the arrangement of the strain gauges com- 





Pt ia 


Fig. 4.3 Calibration of the Dynamometer 
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pensated for any bending or longitudinal straining. Two of the strain 
gauges were in tension and were diametrically opposite while the other 
two gauges were in compression and were also diametrically opposite. 
Bending was compensated for since the gauges were opposite one another. 
Axial stretching was compensated for since two of the gauges were in 
tension and two were in compression. 

Results of the calibration can be seen in Table 4.1 and 
Figure 4.4. The calibration curve is linear and the conversion factor 
from gauge readings to torque, as calculated by the ratio (z MeV; )/(z a): 
is 179 in-1b/mv. 
4.1.3.b Calibration of Twist Measurement 

The lower jaws on the torsion machine were rotated by a 
"dividing head" which was connected by a chain-sprocket mechanism. It 
was found that 5970 revolutions of the handle of the dividing head 
were required to make one complete revolution of the lower jaws. 

The specimen was then gripped in the lower jaws and the spiders 
were attached. The wires from the transducer-lever arm apparatus were 
secured to the outer edges of the spiders and the slack was taken up 
by means of an adjustment in the apparatus. The voltage output readings 
from the linear transducers could then be correlated with true angles 
turned by the upper and lower spiders. Assuming that each revolution 
of the dividing head handle turned the lower jaws through the same 
angle, a calibration curve was plotted which converted voltage output 


from the transducer to angle of twist for each spider. The calibration 






















By 


ntstte sd? to owl .patnisrte [sntbustpnel 1 pntbned uns “ot bagseneq 
xonto sit oftdw adfeogge yi lsolssennth aia bre netenst at “one 
.9t Faoqqo yffsstrtemeth oafs. sow brie notazsquoo nt s19W 2opusp. ¢ 
ordtons sno stteoage ev3v 2epvee edt sont2 WOT batsemeqnes 26w gnre 
nt evew eanusp eft to owl sonte yor beateeneqnos aew pntdoterte fi 
Moheean gino: nt svaw owt bas Hote 

oo sldeT mn? nase 9d "69 nottsrditso oft to etfuee - ay 
vote? noteyvevnes sit bos wesnth ef avi noitetdttso sat bh ; 
aay, T)\C.V,M 3) offset sad yd batsfuolso 25 ,supYOd OF 2pnrbesy splsy 
| ; vita =nt ORE a 
trigmerivensm _t2twT To mots svat i 6) a. a a 


6 yd betstoy stow sf tat ied not2ves ait a0 ews, Vewol saT 


7 
ai z 


2 i 
aL 


$i .metnedoem Jex.01ge-nhedo 6 Yd betdsnnoo 2ew Wofiw “bead enbae ‘a 
:. on 


bsed prtbivib oft to sibnei odd ‘to encitulovey OV@2 tent bavot 
awet wswol sit to noftulovey atefgmeo sno saan od. bor ups 

evebrqe eid brs ewst iswol 982 af Segaisp eng eke nomtoege oT a 
gyaw 2uts sqgs mins yavel-ysoubenstd art mov? 2ettw eAT -bedaests. | 


qu nodst e6w xosl2 sid bas evabiqe aft to 2apbe tedu0 od ot 






zentbssr Sugqtuo spsttoy edT .evtevsqgs ont nt rete ns to 
*lone ound agtw eds * oe bie ) sani f erit me 


eee Pil | be 4 7 saa e 4 “Ty re fi ioe a 4 iti ’ 


a 








a 


79 


readings for the circular specimen and the Square specimen are shown 
in Tables 4.2 and 4.3 respectively. The corresponding curves are 
shown in Figures 4.5 to 4.8. 

There was some slack in the system prior to the calibrations. 
This shows up in the best-fit curves shown in Figures 4.5 to 4.8. It 
took approximately 1/8 of a revolution on the handle of the dividing 
head to take up the slack. 


A curve of the type 


where N is the number of revolutions of the dividing head and V is the 
voltage reading in volts was fitted to the calibration data so that 


the quantity x (N. - cv. ~ a) is minimized with respect to c. There- 
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The angle of twist is then given by 
re 2a ] 
8 = ag7q (cl + 3) ; (4.1) 


Equation (4.1) for the circular specimen became 


6 = 9.15 x 1072 V + 0.132 x 107? rad. 
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for the upper spider and 
2 -3 -3 
8 s58eh8nxni00) eVier0ei32ixelkOencrad? 


for the lower. The corresponding equations for the square specimen 


are 
. -3 -3 
6 2e9s19exthOtheVete0e] 320x010 rad. 


for the upper spider and 


SF LpnOS132 Ixed@ae nad: 


Oreasi22sxwhOs 
for the lower spider. Since the difference of the upper and lower 
angles of twist is required, the term 0.132 x 1073 rad. has no effect 
on the results. 

It was observed that when the transducer cores moved approxi- 
mately one-tenth of an inch, the Output voltage increased by about 0.6 
volts. Voltage readings were taken on the one volt scale and the 
voltages ran from 0.200 volts to 0.800 volts. When the voltage reached 
0.800 volts the transducer core was reset to its initial position by 
an adjustment in the lever arm apparatus, thus in effect resetting the 
voltage reading to 0.200 volts. By using this procedure more sensitive 
voltage increments could be read from the voltmeters. Another reason 


for this procedure was that the final displacements were larger than 
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the allowable displacements of the transducer cores. 
4.1.3.c Experimental Procedure 

As soon as each particular specimen was calibrated the upper 
jaws of the torsion machine were tightened and the experiments were 
performed. 

A constant tension was maintained in the wires by attaching 
some "dead weights" to the cores of the transducers by means of wires 
and pulleys. The tension in the wires caused a negative moment of 
magnitude 1.79 in.-1b. in the specimen. 

It was also noticed that the angle measuring apparatus was 
not reversible due to friction and possibly some slack in the gears. 
For this reason, readings were taken for loading only. 

Pins were inserted in each end of the specimens to help 
ensure that there would be no slip in the jaws of the torsion machine. 
If there was slip however, it had no bearing on the results because 
the moment measured would still be the true moment and the difference 
in the readings given by the two spiders would still give angle of 


twist per unit length. 


4.1.4 Results of the Torsion Tests 

The output voltage readings for the circular cross-section 
and the square cross-section are listed in Tables 4.4 and 4.6 re- 
spectively. These readings are then converted to moments and corres- 
ponding angles of twist per unit length by the calibration equations 


shown in Sections 4.1.3.a and 4.1.3.b. The results are shown in Tables 
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4.5 and 4.7. The angles of twist per unit length for Tables 4.5 and 
4.7 are calculated by dividing the difference in angles turned by the 
upper and lower spiders by the respective gauge lengths. As was 
mentioned in Section 4.1.2 the gauge length for the circular specimen 
is 1.95 in. whereas the square specimen has a gauge length of 1.75 in. 

Figure 4.9 shows a plot of moment M versus angle of twist 
per unit length 6 for the specimen of circular cross-section. 

The classical elastic theory for torsion of a circular bar 


of radius "a" predicts that 


The slope of the M Vs 6 curve in the elastic region then is 


O|= 


ce 
tran | 


Since the radius of the circular specimen is 0.181 in., the initial 


slope of the curve in Figure 4.9 should be 


n(0.181)" 


3 
2 °) 


(32760x 105)ir=s6"34axql0ie: 


This slope is in good agreement with the experimental data. In calcu- 


lating the dimensionless moments M/2ua®, which are shown in Tables 4.8 


6 


and 4.9, the shear modulus up is taken as 3.76 x 10 psi [20]. The re- 
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lation between M/2ua> and 6a is shown graphically for the circular 
cross-section in Figure 4.10 and for the square cross-section in 
Figure 4.11. The initial slope for the circular specimen is approxi- 
mately 0.825 whereas for the square specimen the initial slope is 
approximately 1.15. The classical theory predicts an elastic slope 
of 7/4 or 0.785 for the circular specimen and 1.125 for the square 
specimen [19]. 

The angles of twist were corrected so that the curves in 
Figures 4.10 and 4.11 extrapolated through the origin. For the cir- 
cular cross-section, 0.02 x 105 was subtracted from the angles of 


3 


twist (6a) while 0.03 x 10 ~ was added for the square bar. 


4.1.5 Comparison Test 


The purpose of the next test was to compare the angle of 
twist per unit length obtained by the spider apparatus with that deter- 
mined by strain gauges. Since strain gauges can be used to find angle 
of twist per unit length for circular specimens, a commercially pure 
Aluminum specimen was machined and annealed as described in Section 
Sa eae 

The test diameter of the circular specimen was 0.371 in. 
and the gauge length was 1.942 in. Four strain gauges were attached 
at approximately the center of the test section of the specimen. The 
strain gauges were mounted at 45° to the axis of the specimen so that 
two of the gauges acted in tension and two in compression when the 


specimen was twisted. The gauges were located at 90° to each other 
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such that each similar pair were diametrically opposite. This arrange- 
ment compensated for any bending or axial stretching. 

The calibration curves for angle of twist determined by the 
Spider apparatus are shown in Figures 4.12 and 4.13 for the lower and 
upper spiders respectively. The calibration readings are tabulated 
in Table 4.10. The equations for converting voltage readings to angles 
of twist are 


3 


O-a8s10.” V4 0.526x 10” rad. 


jap) 
iH] 


and 


é 


CAB ae AVE S26°00 105° rad. 


@D 
iH] 


for the upper and lower spiders respectively. These equations were 
obtained by a least squares approximation. 

Table 4.11 shows the strain gauge readings and the spider 
voltage readings that were taken during the test. The strain gauge 
readings are converted to angles of twist per unit length in Table 
4.11 and the spider voltage readings are converted to their corres- 
ponding angles of twist per unit length in Table 4.12. Figure 4.14 


compares the two graphically. 
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4.2 Tension Test 
A commercially pure Aluminum tension specimen with a circular 
cross-section was prepared as in Section 4.1.2. Figure 4.15 shows the 


specifications. 


Ge75u 

leNvae dia: 
Ue 

0.80" dia. 
t75? 


0.367" dia. 





Fig. 4.15 The Tensile Test Specimen 
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Four strain gauges were attached longitudinally at the center 


of the specimen and were located at diametrically opposite positions 


as in Figure 4.16. 


Fig. 4.16 Cross-Sectional View of Tensile Specimen Showing 


the Location of the Strain Gauges 


The load and corresponding average strain readings are listed 
in Table 4.13. Since the diameter of the test section was 0.367 in., 
the cross-sectional area was 0.106 in.?, Also since [20] v = 0.33 and 


Ea. 10° x 10° psi for commercially pure Aluminum, the deviatoric strain 


can be found from 


* 1-2 ra -8 


where 71 is given in units of psi. The corresponding deviatoric 
stresses and strains are shown in Table 4.14 and plotted in Figure 


4.17. In finding the non-dimensional deviatoric stress, u was taken 


to be 3.76 x 10° psi [20]. 
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TABLE 4.1 
Torque Calibration 


Moment Voltage 
(in. 1b) (mv) 
0 0 
2.20 0-073 
6.61 0.038 
1.02 0.062 
15.43 0.087 
19.84 0.111 
24.25 02735 
28.66 0.160 
oa507 0.184 
37 .48 0.208 
41.89 0.232 


46.30 0.256 
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TABLE 4.2 


Calibration of Angle of Twist for Circular Specimen 


No. of Revs. of Voltage Reading Voltage Increase 
Dividing Head (volts) (volts ) 
Upper Lower Upper Lower 


Spider Spider Spider Spider 


0 0.200 0.200 0.000 0.000 
] 0.287 0.312 0.087 Oe bie 
2 0.394 0.432 0.194 0.232 
a 0.525 0.565 Q). 325 0.365 
4 0.620 0.689 0.420 0.489 
5 Q,. 725 0.217 0. Sk 0.617 
6 0.271 0.351 0.671 0. 75] 
j/ 0.377 0.473 On Liz, 0.873 
8 0.509 0.610 0.909 1.010 
g 0.600 0.732 1.000 1.132 
10 0.753 0.270 bus 5.240 
11 0.260 0.402 1.260 1.402 
12 0.37) 0.53] 5.3 ¥.931 
13 0.480 0.658 1.480 1.658 
14 0.597 6.795 ¥.59/ ¥..795 
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TABLE 4.3 


Calibration of Angle of Twist for Square Specimen 


No. of Revs. of 





Dividing Head (volts) (volts) 
Upper Lower Upper Lower 
Spider Spider Spider Spider 
0 0.200 0.200 0.000 0.000 
] 0.310 QO. 328 0.110 0.128 
2 0.410 0.439 0.210 0.239 
3 0.520 0.562 0.320 0.362 
4 0.635 0.690 0.435 0.490 
5 0.780 0.220 0.580 0.620 
6 0.268 Qrse2 0.668 Ga/32 
7 O35 0.460 Oe /go 0.860 
8 0.510 0.600 0.910 1.000 
9 0.620 0.430 1.020 1.130 
10 0.750 0.268 1-160 1.268 
Ai 0.250 0.400 leca0 1.400 
+2 0.350 0.520 1.350 Teo2a 
13 0.465 0.659 1.465 1.659 
14 0.580 0.780 1.580 1.780 
15 0.700 0.300 1.700 1.900 
16 0.230 0.440 1.830 2.040 
17 0.329 0.566 1.929 2.166 
18 0.440 0.690 2.040 2.290 
19 0.560 0.220 2.160 2.420 
20 0.680 0.348 2.280 2.548 
21 0.800 0.464 2.400 2.664 
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TABLE 4.4 
Experimental Data for Torsion of Circular Bar 


No. of Revs. of Upper Lower Moment Output 
Dividing Head Spider Spider (mv) 
(volts) (volts) 

0 0.200 0.200 0.000 
2 0.205 0.328 0.014 
4 0.241 Ow5l2 0.039 
6 0.267 0.701 0.061 
8 0.294 0.298 0.084 
10 0.314 0.482 0.104 
12 0.329 0.684 0.119 
14 0.340 0.300 0.129 
16 0.345 Or522 G2135 
18 0.350 0.750 0.140 
20 0.360 0.406 0.145 
22 0.360 0.625 0.149 
24 0.368 0.290 OI Ky 
26 0.368 0.530 O=153 
28 0.368 0.748 0.158 
30 0.369 0.399 0.160 
32 0.369 0.639 0.161 
34 0.370 0.280 Oni63 
36 0.379 0.540 0.165 
38 0.380 0.790 0.168 
40 0.389 0.440 0.169 
42 0.400 0.656 02171 
44 0.402 O33! O22 
46 0.402 0.582 0.174 
48 0.417 Orcs OF177 
50 0.420 0.451 0.179 
sy 0.420 0.703 0.180 
54 0.431 0.320 0.181 
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TABLE 4.5 


Results for Circular Cross-Section 


Moment Angle of Twist Difference Angle of Twist per 
(in. 1b) Upper Spider Lower Spider in angle unit length 


(x10? rad) (xt0*? rad) — (x10*? vad) (x10*3 rad/in) 
2, SE ea ee ae 


0.00 0.000 0.000 0.000 0.00 
ooo 0.046 1.047 1.00 0.51 
6.98 0.375 roles SY Zod Vala 
10.92 0.613 4.098 3.49 1g 
15.0 0.860 my at 4.85 2.49 
18.6 1.04 Aso Om Sills 
a jeg 8.87 7.09 3.94 
23.1 Pece 10.6 Ges 4.77 
24.2 Teo3 12.4 Lil sy oe 
Zoe| Lacy 14.3 Weed) 6.62 
26.0 1.46 16.4 14.9 7.64 
26.7 1.46 Igs2 165% 8.56 
ole 1.54 20.4 1839 920d 
27.4 1.54 Fay de 20.8 Oia, 
28.3 1.54 24.1 22.0 Li6 
28.6 55 cOce 24.7 i det | 
28.8 Teo 2-1 2025 HS bes 
2972 256 30.) 2Ge0 14.6 
Gao 1.64 Occ 30.6 152d 
503] 1.65 34.3 3250 10nd 
30.3 aye: 30.0 34.6 UG peo | 
30.6 1.83 so 3055 18.6 
30.8 F285 40.3 38.4 Lo, 
Stel 1.85 42.4 40.5 20.8 
a7 deo9 44.3 42.3 Pa Poot 
JeeU eeu 46.2 44.2 nlad 
pave 2.01 48.3 46.3 py! 
32.4 fol 50.1 48.0 24.6 
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TABLE 4.6 


Experimental Data for Torsion of Square Bar 


No. of Revs. of Upper Lower Moment Output 
Dividing Head Spider Spider (mv) 
(volts) (volts) 

0 0.200 0.200 0.011 
2 0.220 0.330 0.041 
4 0.240 0.474 0.074 
6 0.265 .Ga5 0.110 
8 0.305 0.207 0.143 
10 0.320 0.409 Osl62 
12 ge 0.620 0.176 
14 (- 353 0.259 0.188 
16 0.369 0.478 0.194 
18 0.380 Ay Aes 0.200 
20 0.400 O52 0.207 
22 0.411 0.601 OF21G 
24 0.435 0.230 0.220 
26 0.440 0.473 0.220 
28 0.441 0.721 Q2223 
30 0.468 0.350 0.230 
32 0.468 0.600 0:23] 
34 0.488 0.220 0.240 
36 0.490 0.460 0.240 
38 0.492 0.710 0.241 
40 0.516 0.300 0.247 
42 O.517 0.550 0.248 
44 0.518 0.800 0.250 
46 0.538 0.390 0.252 
48 Q.237 0.650 07253 
50 0.541 0.278 0.259 
52 0.543 0.519 0.260 
54 0.544 0.748 0.260 
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Angle of Twist pitt : 
Upper Spider Lower Spider bats ieecs 1" 
(x10t? rad) — (x10*? rad) 

0.000 0.000 
0.184 L- 02 
0.368 2.20 
0.597 mays) 
0.965 4.99 

iiggdl Au 6.65 

Te22 8.38 

1.41 10.0 

7.55 T2t] 

1265 ac 

1.84 16.0 

1.94 18.1 
Marsh ile) 20.0 
ea 22.0 
Zee 24.0 
2.46 25.9 
2.46 27.9 
2.65 29.8 
2.67 egy 
2.68 33.8 
2.90 35.3 
2.91 Pag | 
2.92 35.59 
Seti 41.0 
Oe LN ia? 
Jaro 45.0 
3.15 47.0 
306 48.9 


Results for Square Cross-Section 


TABLE 4.7 


Angle 


(xi0t? rad) 
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Angle of Twist 
per unit Length 


(x10*? yad/in) 
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TABLE 4.8 


Summary of Results for Circular Bar 





. e M = 
DeiTasnoueade 3 (x10*3) ga(x10*?) 6a-0.02 x 10 3 


2ua (x1 ot) 

0 0.000 0.000 0.000 
2 0.056 0.093 0073 
S O.45¥ 0.202 0.182 
6 0.245 0.324 0.304 
8 0.337 0.450 0.430 
10 0.417 Q2575 0/553 
12 0.478 0.714 0.694 
14 0.518 0.863 0.843 
16 0.542 1.03 1.01 
18 0.562 1220 1:18 
20 0.582 1.38 1236 
22 0.598 T2955 Was 
24 0.610 Te #5 Tvs 
26 0.614 Ts93 19 | 
28 0.634 ea 2.08 
30 0.642 2a29 Aa) 
32 0.646 2.46 2.44 
34 0.654 2365 2.63 
36 0.662 2.84 2.82 
38 0.674 3203 3201 
40 0.678 S221 3v19 
42 0.686 oe! 3235 
aS 0.690 3.56 3.54 
46 0.698 3.76 3.74 
48 0.711 3«93 3391 
50 0.719 4.10 4.08 
52 0.723 4.30 4.28 
54 0.727 4.46 4.44 
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TABLE 4.9 


Summary of Results for Square Bar 


Dividing Head”) oy (xT0") —ea(x10t3) eat0.03 x 1073 





2ua (x10*3) 
—_______ ier wer nner prey 
0 0.039 0.00 0.03 
2 0.146 0.096 0.126 
4 0.263 0.203 0.233 
6 0.391 0-321 05351 
8 0.508 0.433 0.463 
10 0.576 0.598 0.628 
12 0.625 0.771 0.801 
14 0.668 0.926 0.956 
16 0.689 lhalg 1.16 
18 0.711 134 133/ 
20 0.736 1E53 1.56 
22 0.746 1.74 1.77 
24 0.782 1.92 1395 
26 0.782 22\3 2.16 
28 0.793 2.35 2.38 
30 @A817 2252 2.55 
32 0.821 2.74 Zo, 
34 0.853 2e02 2.95 
36 0.853 Sea Saab 
38 0.856 3505 3200 
40 0.878 3.49 Se52 
42 0.881 Sele oa75 
44 0.888 3.94 3.97 
46 0.896 4.08 4.11 
48 0.899 4.32 4.35 
50 0.920 4.51 4.54 
we 0.924 4.72 4.75 
54 0.924 4.92 4.95 
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TABLE 4.10 


Calibration of Angle of Twist for Comparison Specimen 


No. of Revs. of Voltage Reading Voltage Increase 
Dividing Head (volts) (volts) 
Upper Lower Upper Lower 


Spider Spider Spider Spider 
ee ee aes 


0 0.200 0.200 0.000 0.000 
] 0.277 0.282 0.077 0.082 
2 0.379 0.382 0-179 0.182 
2 0.490 0.500 0.290 0.300 
é 0.606 0.630 0.406 0.430 
5 Os aly 0.740 OO by, 0.540 
6 0.230 0.265 0.630 0.665 
7 0.356 05395 02756 Oas25 
8 0.460 0.510 0.860 0.910 
9 0.570 0.632 970 1052 
10 0.679 0.769 1.079 1.169 
1] 0.800 0.290 1.200 1290 
Tz 0.300 0.420 1.300 1.420 
1s 0.420 0.545 1.420 1.545 
14 0.519 0.670 1519 1.670 
15 0.630 Or2le 1.630 1.812 
16 0.744 0.340 1.744 1.940 
17 0.250 0.461 1.850 2.061 
18 0.360 0.590 1.960 26190 
19 0.470 0.723 2.070 a Re, 
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TABLE 4.11 
Experimental Data for Comparison Test 


No. of Revs. of Average Strain Angle of Twist Spider Voltage Readings 


Dividing Head Gauge Reading Per Unit Length (volts) 
(x10*®) (x10*? rad/in) ae Re 
0 0 0 0.200 0.200 
2 43 0.46 0.203 0.280 
4 103 Te ld Os 227 0.443 
6 169 ez 0.250 0.621 
8 231 2.49 0.278 0.203 
10 297 3.20 0.308 0.388 
12 358 3.86 0.350 0.569 
14 422 4.55 0.378 0.769 
16 478 be 15 0.425 0.350 
18 546 5.89 0.460 0.540 
20 607 6.54 0.486 0.744 
22 672 7.24 0.537 0.346 
24 7feyi 7.94 0.557 0.540 
26 812 Su/5 0.580 0.754 
28 880 9.49 0.610 0.362 
30 955 10.3 0.620 0.576 
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TABLE 4.12 
Calculation of Angle of Twist from Spider Voltage Readings 
aes Pm ct Angle ope se 
ncrease ngle of Twis ifference per uni 
eaten iecae (volts) (x10*S rad) sin Angle Sein 
Upper Lower Upper Lower (x10*? rad) (x10*? rad/in) 

0 0.000 0.000 0.000 0.000 0.000 0.00 
2 0.003 0.080 0.028 0.678 0.650 0.33 
4 0.027°°0.243 0.252 2.06 Ae 0.93 
6 0.050 0.421 0.467 3.57 S10 1.60 
8 03078™™0.6039 6.7268 9 5.11 4.38 2.26 
10 OF 108**027838 71.01 6.68 507 2.92 
12 O21 505° 0.969 ale40 8.22 6.82 3.51 
14 OF176> -1. 169m 66 9.9] GS. 4.25 
16 O32257 SOM oO mel e45 0235 4.8] 
18 OF 2005. 540 243k] 10.7 sicel 
20 0.286 1.744 2.67 14.8 lZal 6.23 
22 0.337 1.946 3.14 . 16.5 13.4 6.90 
24 W357 2.140 ms3- 33 en Toe 14.8 1362 
26 053807 "2.354 935550082070 16.4 8.44 
28 OPATOr 2.5620 3.530 oly 17.9 9.22 

30 OF420 (2.776 73.92) 923-5 19.6 10.1 
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TABLE 4.13 


Experimental Data for Tensile Test 





ie Average Strain (e 1) a 
: (x10*°) (psi) 
0 0 0 
50 48 472 
100 98 943 
150 149 1420 
200 203 1890 
250 264 2360 
300 361 2830 
350 932 3300 
400 918 3770 
450 1822 4250 
500 3335 4720 
550 5349 5190 
600 rivals: 5660 
650 10130 6130 
700 12060 6600 
750 16610 7080 
800 20860 7550 
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TABLE 4.14 


Summary of Results for Tensile Test 


al -8 
T= he € = e117 Ges Sa) O11 
(x10*®) (x10*®) 
oie Sk. ee ee 

0 0 
41.8 42.7 
83.6 87.3 

125 133 
167 182 
209 237 
25] 329 
293 495 
335 875 
377 1770 
419 3280 
461 5290 
502 7650 
544 10100 
586 13000 
628 16500 
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CHAPTER V 
SOLUTION OF SOME WORK-HARDENING PROBLEMS 
AND COMPARISON TO EXPERIMENTAL RESULTS 


5.1 Determination of the Work-Hardening Function 


The hypoelastic constitutive equation 


Ding 

JSS 7), 5. | 
pe dig + ANT, (5-1) 
is applied, in this chapter, to Engineering problems involving work- 
hardening materials. The solutions are compared to the experimental 
results obtained for the tension and torsion tests involving com- 
mercially pure Aluminum. 


In applying equation (5.1) to the simple tension problem, 


the resulting equation [14] is 





dine Otere 

der tledrandl 
where T = $71 /eu and e« is the deviatoric strain. Therefore 

dT 
h de . 
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deviatoric stress-deviatoric strain data shown graphically in Figure 
4.17, then h can be obtained as a function of J. However in the 

above form, h is indeterminate at the Origin, and near the origin 

both the numerator and denominator are extremely small resulting 

in a large possible error in h. The following work-hardening function 


does not have this defect; 


The hypoelastic constitutive equation (5.1) for a work- 


hardening material will then be written as 


Ope ; 
(5.2) 


— a = gt [ites 
De 44g * NG) Ts, 
where 0 > H > -1 and the second term on the right hand side vanishes 
when J = 0. 

The slope de/dT is approximated at each of the data points 
in Figure 4.17 by a fourth order Lagrange interpolation formula. The 


values are listed in Table 5.1 and plotted in Figure 5.1. 


Let 


and let a and b take values so that the initial and final g(T) is zero. 










rl . he 
, af 

sip ti nt vilgotdaeye nwode etsb ntevt2 oModsiveb-eeerte oitosstyebe 
Stupid Ar VVIBOTAGE TR Mwone Boh S132 DIMDISTVSO-2ests J vso 


7 
. . 
\ i t , = om ° * A 7 
sit nt vaveawok .b to ortonut & 26 bantstdo sd aso nony ai 


s 
: ‘ - , a ok - 
rtorvyo adt yssn bas . rerio snc. 76 taanioyvetebat ef a gmrrot Svods 


pattivesy fleme yvismayixe s cotentmonab bas “oteromun ont ito { 


rt90U7 irnsbran-doow paftworlot sit i nft-voxyg. aidfe2zoq apisi 6 nt 
é a 

;tostsb afdt oven. don aon 

e -} 7 


+ Tai balboa 
6 ngjtinw ed i+ [ltw detretsm p 
Ty 
~< | 
= 
{ | 
. ‘ * F* 
ro on) "Mave a) oad pee : 
Cc} (—<}H + .*b = 3 
e4 ris 
. ee a pe i aa ee ee a yo r Pa j 
ean2tnay sbt2 boned tnpty sad no mist baosee*siy is “f= < Fh < 
pi ~~ la 


" - . a ae © ee ee . 4 4 ae ee ae. ae © torn f : 
cull hab en: 10 Hwe 76 DOSED qqgs of 7b\ab ogo ted oiT ' 
y ‘) } . 
8 


i x 7 ad io a fj 
oiT .6fumyot noftsloqretar spnsipsd Wbio Minot 6 Ya Vt.) swupt? at 






.[.0 svupra af bessolq brs, [.¢ oTdeT at bagerl abe 28 fev 
t ; 5 ; ‘ - 

4) 4 

7 tel 


ine 


1 ver 
' iyi 
. 7 


116 


A Fourier sine series,” 


15 
knT 
g(T) = ) By, sin(=—) , 
k=1 le 
is fitted to the function g(T). Then 
15 
dew fel kal 
Pitan eee cho eh ) Ae (523) 
aon 
Since J= m 1a 
1 Apa? kad 
Tat aaa ay eee cies in(—) (5.4) 
1+H 3 = 
where 
a = 1.000 
b = 168.8 x 10° 
By = 40.57 
B. Se 1S 
By = 1.577 
Ba = -5,783 
Be =73.170 
Be = 0.3168 
Bo = -0.9106 
Be = 0.1436 
By = 0.1017 


*T is the final experimental value of T and in this particular case 


T, = 670.1 x 107° . 
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B,, = -0.0474 


10 
By 0.2405 
Bio = -0.3466 
B13 = sey! 
Bia ===0 0012 
Bis = -0.0329 


Figure 6.1 shows a plot of the work-hardening function H as a function 


of J. 
Recall that for simple tension 
dias 
pray Tete, 
Therefore 
- 
re wi hseS ON 
cle eb uf kT) _ 
E = | T+Ho: al + ie ars oy (—) ee ) We (5.5) 
1@) 


The stress-strain relationship (5.5) is shown graphically in Figure 


5.2 and compared to the experimental results. 


5.2 Simple Shear 


Consider a material that obeys the constitutive relation 


vit ) 
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For simple shear it follows that 





where 


(=p) 
fa) 
(=) 


2452. The material derivative has been used to arrive at 


and J=T 
equations (5.6) and (5.7). It follows from equation (5.6) and the 
initial condition that T = 0 when c = 0, that T = 0 always. Therefore 
equation (5.7) becomes 

dS 


dale 


It follows from equation (5.4) that 
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The relation between S and c is shown graphically in Figure 5.3. 


5.3 Torsion of a Circular Bar 
In the case of torsion of a circular bar equation (5.2) re- 


duces to 


oT 


sae 5 (1 + H) (5.8) 


where T = S95/2u and 6 is the angle of twist per unit length. Also 


2 


J=T™. Equation (5.8) can be integrated to yield 


t 





“p = al + =. L (4) B, [eos (tt 2 cL ott ie 

Now since the moment is defined as 

a 

M = | nr’ Sd 
e) 
it follows that 
l 2 
Hear f 1) ah . (5.9) 
a a 
2ua : 


Let To represent the shear stress when r = a. Then 
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bie 2 Ta Te knTp 
as aly + 0 1 ) (Z) By [cos ( T; ) = Che . (5.10) 





Also for some particular angle of twist (6a) 








Loser | 
2 Dee is kaT 
(F) = (A)[fat +2 1% + JY (4) BLcos (kth - ig 
a 6a 3 pay KT as Te 
From equation (5.9) it follows that 
Ms 
eel: 2 
= - — | Tat +21? + j (2) Loos (==) = WARE 
dua (6a) A V3 k=] f 


(521) 
15 
- [a + bf - ) Bg sin(<) ale 
k=] f 


For each To: (6a) can be found from equation (5.10). Then the integral 
in equation (5.11) is evaluated by Simpson's Rule with Ty divided into 
10 equal segments. As Tg varies from 0 to Te a corresponding (6a) and 
M/2ua® can be found. Some of these values are listed in Table 5.2 

and a plot of moment versus angle of twist is shown in Figure 5.4 and 


compared with experimental torsion results. 


5.4 Torsion of a Square Bar 


The problem of torsion of a material that obeys the constitutive 


equation 
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where a = aye /k@ was solved in Section 3.3.c and 3.3.d. In this 


section the problem is solved for a square bar where 


The solution proceeds as in Section 3.3.c up to equations (3.15) and 
(3.16) except that - of 


is replaced by (H/2J). 
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Equations (3.15) and 
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Differentiating equation (5.13) with respect to Xo yields 
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equation (5.16) becomes 
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where A=8 
B=y 
C=n 
- By x 
D = + 
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Since equation (5.17) is identical to equation (3.22) the rest of the 
solution is the same as in Section 3 except that the coefficients A, B, 
C, D and E depend on the work-hardening function H which is given by 
equation (5.4). The results are given in Figure 5.5 and compared with 


the experimental results. 


5.5 An Elastic-Plastic Constitutive Equation for a Work-Hardening 


Material 
From the classical theory of plasticity, the flow rule for 


a work-hardening material with a regular yield surface is 
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during plastic flow. It follows that 
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where Ty = 55 5/4u and J = > Teele: The total deviatoric strain in- 
crement des 5 is then given by 
Tae 
dew wend ae ae (5.19) 





¥ ij 2d iJ 


In the dynamically correct form, equation (5.19) becomes 





in the elastic range. 
Consider a material that begins to yield immediately upon 
initial loading. (Commercially pure aluminum is an example). Then 


the constitutive equation 
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ar The dhe 
edt ee ke (5.20) 


applies for the complete range of loading. Now since 


Teodke = Tkedke * Thedke » 
let 
Tate = CH) Typdee = > HM! 
where Se ear JO Pe 
Equation (5.20) becomes 
ah = d!. + H (Mo) Tt; (5.21) 
Dt ij CL : 


where (-H) is the ratio of plastic work to the total work done. 
Equation (5.21) for the elastic-plastic flow of a work- 
hardening material is identical to equation (5.2) for a hypoelastic 


body. 
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TABLE 5.1 


Estimated Slopes of Tensile Results 


T(x10*®) & 
0 1.000 (*) 
41.88 1.051 
83.76 1.075 
125.6 1.121 
167.5 ibid 
209.4 1.620 
251.3 2.728 
293.2 5.641 
335.0 14.42 
376.9 28.74 
418.8 42.55 
460.7 52.99 
502.6 57.23 
544.4 62.86 
586.3 vis aM 
628.2 92.23 
670.1 114.1 


(*) & was set equal to 1.000 at the origin 
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TABLE 5.2 


Results for Torsion of Circular Bar 


Ty(x10*®) —_-ga(xt0*) — (x1o*) 
2ua 
0 0 0 
17.41 0.03563 0.02746 
34.82 0.07220 0.05493 
52.23 0.1087 0.08225 
69.64 0.1455 0.1097 
87.05 0.1834 0.1377 
104.5 0.2222 0.1658 
121.9 0.2612 0.1935 
139.3 0.3005 0.2213 
T5667 0.3421 0.2503 
174.1 0.3896 0.2819 
191.5 0.4473 0.3169 
208.9 0.5204 0.3554 
226.3 0.6158 0.3975 
243.7 0.7460 0.4441 
261.1 0.9375 0.4967 
278.5 1.238 0.5532 
296.0 1.708 0.6053 
313.4 2.401 0.6459 
330.8 3.339 0.6766 
348.2 4.521 0.7046 
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CHAPTER VI 
CONCLUSION 


6.1 Estimation of Errors 
6.1.1 Errors in Numerical Procedure 
For both the elastic-perfectly plastic and the work-hardening 


body, the equation 














2,3 2/3 23 : : 
; 9 (=e) 3 (3) 3 (<5) (a) a(s5) a) ae 
ax,” 3X1 9X5 aX5° Oxy te 


was solved for (38) at each angle of twist per unit length in the 
torsion problem. A new stress function was calculated each time by 
applying the formula 
p(erae) = y(6-Ae) + 2A0(S%) (6.2) 

at each point in the square cross-section. A corresponding moment 
was then found by integrating the stress function over the cross- 
section by using Simpson's Rule. 

An estimate of the truncation error introduced when replacing 
equation (6.1) by a finite difference approximation is obtained by 


changing the mesh size. In Chapter III and Chapter V a mesh size of 


1/10 is used. The problem was also solved with a mesh size of 1/6. 
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It was found that a local truncation error of 0.2% is present in the 
calculated value of the rate of increase of the stress function at the 
center of the cross-section. After 800 increments in angle of twist 
per unit length this truncation error had accumulated to a value of 
0.7%. These truncation errors are at least 1000 times as large as 

the error involved in using Simpson's Rule to find moments from the 
stress function. 

The truncation errors in the stress function, as described 
above, caused a local error in the moment of 0.13%. After 800 incre- 
ments in angle of twist the error became 0.3%. 

The increment in angle of twist per unit length was then 
doubled for the case of the mesh size equal to 1/10. This made no 
Significant change in any of the results. The computer program 
printed out results to eight significant figures and after one thousand 
increments in angle no change could be noticed in the eighth digit of 
the stress function or the moment. 

Roundoff errors were minimized by using double precision 
accuracy when making the computer calculations. A computer library 
subroutine (DGELB), which is designed specifically for inverting 
matrices with a banded structure, was used to solve the 55 equations 
at each angle of twist per unit length. The solution of these equations 
was carried out with a precision of 107 !8 [244]. 

In addition to the above, there are errors in the calculated 
moments due to errors in the work-hardening function H. 


H is found from the slope of the stress-strain data from the 
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tension test as described in Chapter V. From equations (5.3) and 


(5.4) 
doe al ree Pia 
dr > isn 74) 8 m vee i Bk sin(———) . (603) 
k=1 We 
Consider 
d 
(se) = a le e,[b/ = J - J Bk sin(K™) J (6.4) 
1 1 Je 
and 





15 
de ] 4 sr Kivd 
(==) = = | *+>e [by = J - Bk sin(——) ] (6.5) 


iP 
where 0 < e, < 1 and ey > 1. Therefore 
de de de 
(Gt), Sat < (ar), (6.6) 


where e are the calculated values from equation (6.3). It follows 


then that 


9 6.7) 


where the equalities in equations (6.6) and (6.7) apply only when 


J = 0. During loading, for a positive angle of twist, T13 Ztu 
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There- 


it follows that Mo <M, where M is the moment and M, and M, are the 


moments calculated when Hy and Ho are substituted respectively for H. 


Equations (6.4) and (6.5) are imposing an error on de/dT - 1 


for all_J and yet the functions Hy and Ho are still bounded by the 


inequalities 
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If the function de/dT - 1 can be determined from the experimental 
tension data to within the imposed error at all points then M, and 
My are certainly upper and lower estimates respectively for the 
moment. 

It was found that M, and My for the problem of torsion 
of a square bar varied only about 2.5% from the previously calcu- 
lated moments M when a 20% error in de/dT - 1 was imposed by letting 


e, = 0.80 and e, = 1.20. 


2 
During the experimental tension tests, assume that the 

strains € were obtained with an accuracy of + 3%. The slope of the 

T vs ¢ data was obtained at each data point by a computer library 

subroutine which used a fourth order Lagrange interpolation formula. 

An estimate of the maximum errors in the calculated slopes de/dT [25] 

is given below. 
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where e. is the error in the ordinate e« and AT is the equally spaced 


intervals in stresses T. 
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15 
de anes ead 
=~ = 1+ bT - } Bk sin(=—) 
gi k=1 "e 
the truncation error is given by 
415 4 
A kt kaT 
Qer =- a, ) (=) Bk sin(=—) 
Ts Peet Ts 
Therefore since Te = 16A , 
4 15 
ep = = OD. Ok Bk sin (6.8) 
30(16)" k= f 
and 
os) 0.03 
Cl ee (6.9) 


nie She 


where ey and Ep are the maximum truncation and roundoff errors. Table 
6.1 lists the slopes de/dT and the corresponding errors calculated at 
each point from equations (6.8) and (6.9). It appears that the error 
in (de/dT - 1) may be approximately 100% when de/dT is small. However 
as the slope increases, a 20% error may be assumed for the purpose of 
establishing bounds on the error in moment. The upper and lower esti- 
mates, M, and Mo » were calculated assuming 100% error in de/dT - 1 when 


<3 xX 1078 and a 20% error when J > 3.x 1078. The factors e, and e, 


Se Ree te Red ee IE a. Dh a ee ee, TO 
*It may be noted that the fifth derivative in equation (6.8) is the 
fifth derivative of a best-fit curve through the calculated slopes. 
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in equations (6.4) and (6.5) take the following values; 


0.10 i) ep yer 
e = 
-8 
0.80 (i> 3 x 100°) 
2.00 (eo 107) 
&5 = x 
1.20 eee Gas) 


The corresponding curves for Hy and Ho are shown in Figure 6.1. Table 
6.2 shows the percentage errors in the moments as calculated by the 
equation 


(M,-M) 
(upper/lower) bound = yw X «(100 (i=1,2) 





for a square bar at various angles of twist. It can be seen that the 
maximum difference between the calculated moments M and the moments M, 
and My is approximately 6% and occurs at a relatively small angle of 
twist. The difference then tends to decrease towards approximately 
2.5% as the angle of twist per unit length increases. It appears 
that the large possible error in H for small J has little effect at 
larger angles of twist. 

The corresponding errors in moments for the circular bar 
behaved in the same manner. A table for the circular bar similar to 


Table 6.2 is not included because when H changes, the angles of twist 
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as well as the moments change and the moments could not be compared 
at any specific angle of twist. 
6.1.2 Experimental Errors 
An estimate of the errors for the tension test is given in 
Section 6.1.1. 
The following discussion on errors for the experimental 
torsion problem applies to both the circular and square cross-section. 
The non-dimensional moment is given by M/2ua>. The moment M 
can be obtained from the expression M = 179 V to within 1% of the true 
moment where VY is the voltage reading in m.v. and M is moment in units 


oe within 1% 


of in.-1b. Since "a" is measured to within + 0.001", a 
of the true value. Also it is assumed that u is correct to within 1%. 
The fact that the initial slope in Figure 4.9 compares so well with the 
classical elastic theory for torsion of a circular bar indicates that 
p= 3.76% 10° psi iS a reasonable value to use. Therefore the ex- 
pression M/2ua> may have a maximum error of 3%. However, the most 
probable error is approximately 1.7%. 

Originally, neither of the experimental curves for moment 
versus angle of twist went through the origin. This can be seen in 
Tables 4.8 and 4.9 where a constant was added to the 6a values. This 
may partially be due to a small moment induced in the specimens when 
the upper jaws of the torsion machine were tightened prior to the experi- 
ments. It may also be caused by one of the spider mechanisms having a 


different amount of slack than the other. 
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Section 4.1.5 compares the angle of twist measured by the 
Spider apparatus to the angle of twist measured by strain gauges. 
Figure 4.14 shows that if a constant is added to the angles of twist 
determined by the spider apparatus, the two different methods give 
almost identical results. 

From the above discussion and the close agreement of the 
calibration data with the calibration curves as shown in Chapter IV 
the author makes the assumption that the angles of twist (6a) are 


accurate to within 3%. 


6.2 Discussion of Errors 

It appears that the numerical calculations of the theoretical 
moments for the square cross-section could be improved by decreasing 
the mesh size. However, the numerical technique seems quite adequate 
until the error in H can be reduced for the case of the work-hardening 
solid. 

The theoretical and experimental torsion results compare 
quite favorably. The theoretical moment curves for both the square 
and circular cross-sections deviate from the experimental results in 
the same manner and at approximately the same angle. The angle at 
which the greatest deviation occurs is approximately the same angle 
at which an error in H would cause the greatest error in moment, as 
seen in Table 6.2. 

The calculations shown in Chapter III for the torsion problem 


of an elastic-perfectly plastic material would not involve errors due 
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to an estimated work-hardening function. 

Error estimates in this chapter are concerned only with the 
torsion problem because this is the only problem to which results could 
be compared. An error analysis on the numerical methods used for the 


other minor problems in this thesis can be done using similar reasoning. 


6.3 Discussion of Results 

In Section 6.1.1 it was shown that if Ho SH for all J, the 
moment My calculated by using Hy in the hypoelastic constitutive 
equation would be a lower estimate for M. The same reasoning can be 
applied to the elastic-perfectly plastic results shown in Chapter III. 
It follows then that hypoelastic theory gives an underestimate for the 
torque of the classical elastic-perfectly plastic solution for the 
torsion problem. This is also true for the Simple shear and the simple 
extension problem. 

In Chapter III the simple shear problem was solved with the 


assumption that the material obeyed the constitutive equation 
Bilerc oe appre (6.10) 


where a = Que /ke and M' = Thode: 

Truesdel1 [13] solved the same problem but interpreted the 
results differently than did the author of this thesis. The constitutive 
equation used by Truesdell [13] is identical to equation (6.10) except 


that the constant a® is replaced by 1/K-, Therefore 
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From the results of Section 3.2 and the comments made at the end of 
that section it is clear that if the dynamically correct form of the 
stress rate is used the solution of the problem will consist of a 
shearing stress S and a tensile stress T. The ratio T/k approaches 

k/u asymptotically from below while S/k rises above J - Ke rue and 

then approaches A ~ Kee asymptotically from above. The point at 
which S reaches a maximum is defined as the hypoelastic yield point 

by Truesde1] [13] and the von Mises yield is obtained when T2 + S@ = k2. 
It is shown in Section 3.2 that the asymptotic values shown above are 
valid for the hypoelastic theory and the classical elastic-plastic 
theory if the von Mises yield criteria and the Jaumann stress rate is 
used. However in the hypoelastic theory (12+ 52)/k* approaches unity 
from below while in the classical elastic-plastic theory (T+ S°)/k2 = ] 
once the yield point is reached. If the material or local stress rate 
is used, both of the above theories predict that T = 0. 

Truesdell [13] compares the results of the hypoelastic theory 
in which the Jaumann stress rate is used to that of the classical 
elastic-plastic theory in which the material rate of stress is used. 
These comparisons lead to various incorrect deductions. 

As is mentioned above, a normal stress T appears in Truesdell's 
solution because the Jaumann stress rate is used rather than the ma- 


terial derivative. 
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It can be shown that the hypoelastic yield always occurs at 
a smaller strain than does the von Mises yield as predicted by elastic- 
plastic theory. Truesdell deduces then that since the hypoelastic 
yield is a maximum for shearing stress S, it is impossible for the 
Shear stress to increase continuously to the von Mises yield as ob- 
served by T.Y. Thomas [22]. Although this statement is true, it is 
trivial because S is merely a measure of shear ona particular plane. 
From the Mohr circle shown in Figure 6.2 it is clear that the normal 


Stress T induces a further shear at each point in the body. 


Shearing 
Stress 


aa 1 






Normal Stress 


Fig. 6.2 Mohr Circle for the Simple Shear Problem 
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The Mohr circle in Figure 6.2 is for a typical element in the body as 


shown in Figure 6.3. 


oF. 





Fig. 6.3 Typical Element in a Body During Simple Shear 


The maximum shear stress t is found from the relationship 


Sle tes 


and it was shown in Section 3.2 that t approaches k asymptotically 
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from below. 

Simple shear consists of pure shear with a rotation. If the 
problem of pure shear was solved, no normal stress T would appear since 
the Jaumann stress rate would reduce to the material rate of stress. 

The resulting shear stress for a shearing angle of /2 would be identical 
to the maximum shear t for the simple shear Problem at a shear angle 

of ¢. Therefore the normal stress T arises in the simple shear problem 
due to a rotation through the angle $/2. 


For most metals ae 10°. 


Therefore K is very small and for 
the problems solved in Chapter III, the stresses get very large be- 
fore the rotations become significant. One is then justified in using 
the material derivative in place of the dynamically correct stress rate. 

Truesdel] [13] also concludes that "since for very small K, 
the shear stress at hypoelastic yield is substantially the same as the 
Shear stress at the von Mises yield, the shear stress as a function of 
Shear strain will rise quickly to a value approximately K/v2, overshoot 
to a slightly larger value, and then level off again". Regardless of 
whether Truesdell is considering the shear stress S or the maximum 
shear stress t, this conclusion is erroneous for any permissible 
value of K. Section 3.2 of this thesis shows that t approaches k from 
below. Therefore t/2u approaches k/2u or K/v2 asymptotically from 
below. 

Chapters IV and V consider the experimental measurement and 


the theoretical prediction respectively of problems involving a work- 
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hardening material. The results compare favorably and it is felt that 
the method described in Chapter V may be used quite confidently to 
predict the relation between M/2ua® and 8a for a work-hardening ma- 
terial with a square cross-section. Since angles of twist are ex- 
tremely difficult to measure for a square bar, it would be a simple 
matter to perform a tension test on the material and feed the resulting 
stress strain data into a computer program which would predict the 
performance of the same material under torsion. This technique, of 
course, is not restricted to the torsion problem. Torsion of a square 
bar was chosen as an example since experimental results are sometimes 
difficult to obtain. The results for simple shear, shown in Figure 


5.4, would hopefully compare quite favorably with experimental results. 
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TABLE 6.1 


Estimation of Maximum Errors in Slopes 


de Truncation Round Off 
dT Error Error 
1.000 0.000 0.000 
1.051 0.00894 0.0460 
1.075 -0.0150 0.0940 
era 0.0292 0.143 
aes -0.0130 0.196 
1.620 -0.00221 ORAS 
2728 0.187 0.354 
5.641 ~0.160 0.533 
14.42 = 375 0.942 
28.74 0.381] Teo 
42.55 “0.015 S )e) 
D2299 0.662 5.69 
STE: -0.0958 8.24 
62.86 -0.638 10.9 
7597] 0.616 14.0 
Uzecs -0.548 1756 
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TABLE 6.2 


Estimation of Maximum Errors in Moments 
a2 Te OT Raxiium Errors in Moments 








Ga (M, -M) (M.-M) 
(x10*3) M x 100 aria x 100 
0.000 0.0 0.0 
0.164 fae -4.0 
0.329 6.0 -5.8 

493 5.4 -5.5 
0.657 4.6 -4.7 
0.822 4.0 -4.0 
0.986 Sef =3.6 
Tis 3.4 =3.3 
Ree) ged -3.] 
1.48 Shy -2.9 


oe 







8.0 3J8AT 


zdosmol nt e10173 muaixstt to makaeaitieag = 





001 ae 00! al ad oe | ie 
ae ‘<1. aia) gal 


15] 


O02 


WwNULWN|Y d4Ndg ALL ePLOVaWWO) UOJ UOLZOUNY BULUapuRH-y4OM [°9 “BL4 


i 


8l 9 vl Cl Ol 


8 


On = 


80- 


90- 





tunimufA atu vilstovermod rot nofsanud oninsbysh-s9Ow Wie 


ref 


152 


BIBLIOGRAPHY 


Truesdell, C., The Simplest Rate Theory of Pure Elasticity, 
Comm. Pure Applied Math., 8; p. 123-132, (1955). 

Prager, W., An Elementary Discussion of Definitions of Stress 
Rate, Quarterly Journal of Applied Math., 18; p. 403, (1960). 

Noll, W., On the Continuity of the Solid and Fluid States, 
Journal of Rational Mechanics and Analysis, 4; p. 3-81, 
(1955). 

Jaumann, G., Geschlossenes System Physikalischer und Chemischer 
Differenzialgesetze, Sitzber. Akad. Wiss. Wien. (Ila) 120; 
p. 385-530, (1911). 

Rivlin, R.S., Further Remarks on the Stress-Deformation Re- 
lations for Isotropic Materials, Journal of Rational 
Mechanics and Analysis, 4; p. 681-702, (1955). 

Bernstein, B., Hypoelasticity and Elasticity, Arch. of Rational 
Mechanics and Analysis, 6; p. 90-104, (1960). 

Bernstein, B., Relations Between Hypoelasticity and Elasticity, 
Trans. Soc. Rheol., 4; p. 23-28, (1960). 

Koiter, W.T., General Theorems for Elastic-Plastic Solids, 
Progress in Solid Mechanics, p. 167, Amsterdam, North 


Holland, (1960). 


gar | mani - 





















YHIARDOI 818 


Nib eera 9709 to yrosdT steh tzatqnre sat ,.) ,fhebasuiT Ff 
(2201) SEF-EST .q 28 «ASS bePlgg® gut samod 

22a9t2. ho en ottintts i to porezvoera \wistnamet 4d nA ..W .Yope1d «=. 7 

(0081) .€08 .q 28 « Adem ba fadA to Tsirivel yinetagup ,otsA - 

.zatete otufy bas brioe ont to wtunttned siren ,.W, fon ce 

 FB-E) «gg ‘ohaute bis 2othsdoeM fenotsea to hemwat 


_ »(22e1) 


‘ | 


NSdoetmed) bau wiloe Feared mateye eanae20 lipase | efinsmusl vB | 
sOSf (sll) .notW .22W .bsxA weedsst2 _ostszeplstxnanettto > eee 
(TIE) 088888 ..g 

-sfi nofisnivotsd-2esite edt no zdvemeh edtyen eee entivist ee 

(snotsef to TsnyucG , 2fstiots staovteel yar anette ' 

(2827) ,800-189. .q ah, etayenh bas eotnarooN 

lanoftsA to-.AovA ,\ttofsesta bas Yiotvesisoqyn | Bf ae i) 
080) , S07 -O€ .9.79., 20 2uh enh bas 2atnediogM 

atoltesls iad 8 roides toca nsewitod nots teh ~-8 oti 


ES Ona. 
sG81) -88 al = toe Candi x 





rn i 
- Ay 
_ ; : se 7 = 
7 or }Oc a ade ae 1 TOT 2 >a" c aT 


bane 


10. 


hh. 


tex 


13. 


14. 


Use 


16. 


Fis 


18. 


19. 


20. 


183 


Hill, R., The Mathematical Theory of Plasticity, Oxford University 
Press, Oxford Engineering Science Series, (1950). 

Thomas, T.Y., Plastic Flow and Fracture in Solids, New York; 
Academic Press, p. 94, (1961). 

Green, A.E. and Adkins, J.E., Large Elastic Deformations and 
Non-linear Continuum Mechanics, Clarendon Press, Oxford 
(1960), p. 340. 

Truesdell, C., Second Order Effects in the Mechanics of Materials, 
International Symposium on Second Order Effects in Elasticity, 
Plasticity and Fluid Dynamics, Haifa, b962)Rip slo. 

Permagon Press, (1964). 

Truesdell, C., Hypoelastic Shear, Journal of Applied Physics, 
27; p. 441, (1956). 

Green, A.E., Hypoelasticity and PlasticitygaProc! IRoy: Soc. % 
A234; p. 46-59, (1956). 

Green, A.E., Hypoelasticity and Plasticity II, Journal of 
Rational Mechanics and Analysis, 5; p. 725-734, (1956). 

Ramberg, W. and Osgood, W.R., Description of Stress-Strain Curves 
by Three Parameters, N.A.C.A. Technical Note No. 902, (1943). 

Saint-Venant, B., Mém. Savants Estrangers, vol. 14; (1855). 

Prandtl, L., Physik. Z., vol. 4; (1903). 

Timoshenko, S. and Goodier, J.N., Theory of Elasticity, McGraw- 
Hill, New York, (1951). 

Metals Handbook, 1948 Edition, The American Society for Metals, 
sgn oa LE 
























bef | | ets ai 


oy 

me: 

a ; 

yitovevtnl brotxO .Vgfott26(9 to yiasdT Isoftamanssm ent ‘a neene fe ay 

.(0d2T) ,e29r4s2 goneioe ontreantend bvoitx0 ,e2aaTd © | 

-170) wall seb} foe nf exude’ bis wor? 279 2679 gs ¥-T  eomonT a 

-(13eT) ee. 6, 228849 stmabsaAeis” 6 | | 

bis enottsmotsd sttzals sp183 Ah “2h bre .3.A . e918 a) 

byotx0 , 22079 nobnaisld., 2orneriogl yun aod meont{-noAsy yf | 

OBE gy (0881) 

_2lsiveicm to 2ornerlosM ont af etosttd yeba0-pneaee gay . i fobeauyt 
NI totseslsaeat at59%t3 iaby0 bnoose no mrf2cqnye feaokienvadn! 
of .q .(Saer) atten < 29 rms ny bruld. bas wrottes lS 

_(naer) 22919 HoPgsiede oy Ma 

,2atevnad beflqod %X renidob seeene sth2eTeogyy ad gi febesurT a | 

(G22T) TRA ove g¥8 | 

~-202 VOR .d079 ,Vsroiseald bas yi tot zsfoequH .-3.A pneerd ar i 

A820) .2eOb gqaMOSA OF 

to fanwol. 11 ytioites(9 bos VM Mes isc AG nseT0 ef 

(OOF) AEX~ASY .q pe pareylenA bas 20 :6il99M snorted) rt 

Zevud nisi32-22exd2 to na ttqiaesd 2 fl W ehoogat bas .W cohenidal: j 

46220) gS0@ .o Sioh (sotndoal. .A.d, h. “ natenere® ast ‘di “7 

raulae ihf i fov meee Posi 5 lage 


‘ 
a iv - els cable sh FP 
a he i oe sath 7 ol ee 7 we, ee, F Ad ty 


om 
» et Ss 


i 


ert 
23. 


24. 


25. 


26. 


er. 


28. 


the 


Esser, Hans, and Ahrend, H., Kann die 0.2% Grenze durch eine 
libereinkommenfreie Dehngrenze ersetzt werden, Arch. f. 
Eisenhiittenw, vol. 13, no. 10, p. 425-428, (1939-40). 

Thomas ,°T:Yz, Proc. Natl. Acad. Sci.,°413°p. 720; ° (1955). 

Saint Venant, B., Comptes Rendus Acad. Sci. Paris, 70, (1870), 
p. 473; Journal Math pures et app., 16, (1871), p. 308; 
Comptes Rendus Acad. Sci. Paris, 74, (1872), p. 1009 and 
1083. 

Subroutine DGELB from I.B.M. System/360 Scientific Subroutine 
Package. 

Subroutine DDET5 from I.B.M. System/360 Scientific Subroutine 
Package. 

Froberg, C.E., Introduction to Numerical Analysis, Second 
Edition, Addison-Wesley Publishing Company, (1969), p. 336. 

Ralston, A., A First Course in Numerical Analysis, McGraw-Hill, 
Inc. (1965). 

Lanczos, C., Applied Analysis, Sir Isaac Pitman and Sons Utd, 
(1957). 

Greenberg, H.J., Dorn, W.S., and Wetherell, E.H., A Comparison 
of Flow and Deformation Theories in Plastic Torsion of a 
Square Cylinder, Plasticity: Proceedings of the Second 
Symposium on Naval Structural Mechanics, Brown University, 


R.I., (1960), p. 279-296. 


154 


bai 


ents fob asdeid £8,0 orb nner 
Webyee Tsdyverg ee ee 

lov wieasirinsette ” 
.¥.T . 2smonT 


+ .ow 


PONSQECT) ~25R-298 aq, Of .on eT 


.(deeT) (OST .d gtk, oe VbeSA UTTeR OeIaT | 


(ON8T) ON ,2tts9 . fot .beoA cubash e9agnio8 

to zevoq Aden fsnuee - 

. beoA ~ Bibi” 2i8.qnI0d 
:  EBOT 


MALL mon? 9.13380 ont tvordue 


Leeks (TBP) Ol ¢. Gas TEXAS oq 


bos @COOT .q .{SV8T) . AN .2misd . roe 


mtipordue oftiinsroe 08b \ns2 2ve 


< | -9psAos4 
oftrinstoe Oat \neteve 


snriuordue 


.M.G.2 movt 87200 antiuordue 
i ae 2p54989 
¢ 64.2 yprodin3 


-\asquold pordetl dud ystesW-no2tbbA ,norsiod: 


broose ,2zfeylenA [saris oF noftouborspl 


966 .q , (R0GT) 


_ Ti tH-wsvooM ,2fevfsoA [satya nt seswed S200 A 
(9 ,20soneJ 


.-btd enoe bose asmtr9 oss2l +62 .2feylanA befiqas 


(SCT }r oe? ¢ 


noetraqito? A. 4.13 ee bas eal rod , 0.8 ern 
pen geri aati t nt eines vot sevahsa ed wory to ine 


ad j a 


pam 4%) wat 


a - 
_~ 


LM Brie bas , etal feet ie 


a 
iy 
€s 


«8 ,tnengV totse 


, A od 2teA 



















8 
ti if 










ae 


155 


30. Herakovich, C.T. and Hodge, P.G. Jr., Elastic-Plastic Torsion 
of Hollow Bars by Quadradic Programming, Int. J. Mech. SC, 
U01969) 5" pt 53-63. 

31. Leigh, D.C., Nonlinear Continuum Mechanics, McGraw-Hill, New 
York, (1968). 

32. Hodge, P.G. Jr., Herakovich, C.T., and Stout, R.B., On Numerical 
Comparisons in Elastic-Plastic Torsion, Journal of Applied 


Mechanics, 35, (1968), p. 454-459. 









eer 


norevoT att2sf4-ottesfa ..1b 13.9 .spbotebae ,7.0 ,dotvedetel 
.«fae .toph .6 tol ,onimmeipey9 otbsibsup Yd 296d woll on to il | 
£9262 .q (000r) Me 
wail | PEE wsrxdoM .2otnetoah muuattiod ssantingd ,.9.d .dpted fe 
| RRA AO ee ts 
[sorvanul. v0 OF pduod2 DNG’ «led slo tVORESH igh 2.9 ,9pboH iss 
baffagA to Isnyuol ,nofetor si}esTGno td 26 f3-at enoetieqnod = i 
P2d-Rah .q .(BOEP) ~eE eahmeroem 


7 





a ae 
* > a. et pera ’ 
: Z ; ve wre a | i 
I xe 4 , Y 

















